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Abstract 

The purpose of this project was an evaluation of the 
pentormancer orsttheatdeadbadtmyDahlininand fSmithyi predictor 
algorithms for bottom composition control of a binary dis- 
tillation column. The control behaviour achieved using these 
algorithms was compared to that obtained using discrete 
proportional-integral (PI) and proportional-integral- 
derivative (PID) alyorithms. 

Although the distillation column exhibits highly non- 
linear dynamic behaviour as well as time delay, first order 
plus time delay transfer function models were used to repre- 
sent the dynamic behaviour of the distillation column for 
use in the design procedure for the deadbeat, Dahlin and 
Smith predictor algorithms. 

Simulation studies of these algorithms were carried out 
uSing a nonlinear distillation column model. The simulation 
results showed that for the case of very small model errors, 
the deadbeat algorithm, the Dahlin algorithm and the Smith 
Puedietonfiwere psuperior to the PI and PID algorithms for 
servo control operations. However, the performance of these 
algorithms, tuned by uSing conventional tuning methods, was 
inferior to that achieved using the PI and PID algorithms 
for the feed flow rate disturbance. Subsequently, = an 
improved tuning procedure was utilized and with this proce- 
dure, the column control performance achieved with the dead- 
beat, Dahlin and Smith predictor algorithms was superior to 


that obtained using the PI and PID algorithms. 
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The experimental evaluation of the different control 
algorithms was conducted on a 22.9 cm diameter, eight-tray 
pilot scale distillation column unit interfaced with a dis- 
tributed system of HP 1000 computers. An on-line process gas 
chromatograph was used to provide direct measurement of the 
bottom composition. The primary investigation was the per- 
formance of the Dahlin and Smith predictor algorithms for a 
-25% step change in feed flow rate to the column. From the 
experimental results, it was found that the control perform- 
ance achieved using these algorithms tuned by using the im- 
proved tuning technique was much better than that of the Pl 
algorithm whose poor performance was due to the significant 
amount of time delay present in the process. Use of the PID 
algorithm, gave performance that was close to that achieved 
using the Dahlin algorithm and the Smith predictor. For ser- 
vO operations, employing the tuning constants established 
for regulatory contro!, the Dahlin algorithm provided the 
Deotmecaidumcneaumost  CONSiStent, ‘performance for +12.) step 
Changes in set point. 

It should be noted from both the simulation and experi- 
mental results, that the significantly improved performance 
which resulted using the PID algorithm compared to that with 
the PI algorithm can be attributed to the presence of the 
derivative action. So for controlling time delay processes, 
Since the PID algorithm is model independent, it should be 
used idre sirequently for industria! applications than at 


present. 
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Tq = 187.0) 


Dahlin Algorithm Control of Bottom Composi- 
ErongehforerlZeSteptChangerin Set “Point 


@Bun? SSDAHO1L Sania. ¥5 24), Tp = 563305 
Tae 182.0, i = PA.) 


Dahlin Algorithm Control of Bottom Composi- 
tion for -1% Step Change in Set Point 

(Run S-DB02; K, = -4.29, Tp = 633.0, 

Var eto? 0k p= e20 00) 


Smith Predictor Control of Bottom Composi- 
tion for +1% Step Change in Set Point 
(Run S-SLPO1; Ky = -5.24, Tp = 563.0, 
We eae ee = 0638, Ky = 07220) 


Smith Predictor Control of Bottom Composi- 
ELOn@fone=lZeStepeChangenin Sat FPodnt 

(Run S-SLPO2; Kn = -4.29, Ty = 633.0, 

ere tor soe ke = 0.945, Ky =h-0.263) 

PI Algorithm Control of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 
CRUM S=PLOS;MK OG £9,245) K7.=7-02064) 


PI Algorithm Control, of Bottom Composition 
for a +25% Step Change in Feed Flow Rate 


(Run S-PI04; Ko = -0.267, Ky = -0.133) 


PID Algorithm Control of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 
(Run S-PID033°KE =°-0:262, Kyp°=*+0.1015 

Ky = -0.295) 


PID Algorithm Control of Bottom Composition 
for a +25% Step Change in Feed Flow Rate 
GRUAASSPLOOS SOK. Contge3igf Keksoeol prt 

c at , 
Kp = -0.324) 


Deadbeat Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Rate (Run S-DB03; K, = -4.29, Tp = 633.0, 


p 
Tq = 187.0) 
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Deadbeat Algorithm Control of Bottom Compo- 
Sition for +25% Step Change in Feed Flow 
Rate (Run S-DBO4; Ky = -5.24, Tp = 563.0, 
ae e870!) 


Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Rate (Run S-DAHO3; K, = -4.29, Tp = 633.0, 
nea wee eek = 20.05 


Dahlin Algorithm Control of Bottom Compo- 
Sition for +25% Step Change in Feed Flow 


Rate (Run S-DAHO4; K, = -5.24, Tp m5 05.0), 
Ten 7.0 s An = 28.0) 


Smith Predictor Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Raemch Utiges 5. P00 sek weet = 4229, Ty = 633.0, 
tid od AGS MOA al pas 0.856, Kr = -0.911) 


Smith Predictor Control of Bottom Compo- 
Si tion) for) +25% Step Change in Feed Flow 
Ree emenunes so LPOds Koes 6 24. Tp = 563.0, 
Lied Cehlncn ET DS tell wp SRE Rr 0639 


Deadbeat Algorithm Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 


Rate (Run S-DBO5; Ky Soe ighte, Wns 456.0, 
TPF 877°50') 


Deadbeat Algorithm Control of Bottom Compo- 
sition for +25% Step Change in Feed Flow 
Rate (Run S-DBO06; Ky Sail 8, iy S Sees 


iy = 87720) 


Dahlin Algorithm Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Rate (Run S-DAHO5; K, = -3.12, Tp = 456.0, 
Te =? 1987740 (2 C= 167.5) 


Dahlin Algorithm Control of Bottom Compo- 
Sition for +25% Step Change in Feed Flow 


Rate (Run S-DAHO6; K_ = -1.87, Ty = 352.0, 
ear) 3 P10. eye eos, OF) 


Smith Predictor Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Rate (Run S-SLPO5; BS Hliee er aso. 0, 


Ky = 
Nic eee Sb a AUR LS -0.856, K; = -0.911) 


Smith Predictor Control of Bottom Compo- 
sition for +25% Step Change in Feed Flow 
Rate (Run S-SLP06; Kp Se ovem (p= 352 60, 
esto, Kees -0.547, Kris -0-623) 
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Effect of the value of A in the Dahlin 
Algorithm on the Response of the Controlled 
and Manipulated Variables for a -1% Step 
Change in Set Point 


Effect of Algorithm Constants in the Dead- 
beat Algorithm on the Response of the 
Controlled and Manipulated Variables for a 
-1% Step Change in Set Point 


Schematic Diagram of the Pilot Scale Dis- 
tid katikon Co umn Win t 


PI Algorithm Control of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 
GRUMy B—PS24352 KK l= &0:56.014,, Kiy, =* =0.2095,) 


PobwALoOrdc nim COnurowsor BOLtOm Composition 
for a -25% Step Change in Feed Flow Rate 
(Rum ¢2 <PIDD 2itesk Cherator. 523 stk, =) 120.117, 

Kp = -0.200) 


Dahlin Algorithm Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Rate (Run E-DAH13; CoN 3623, WAP, al p-t=, 19.18,-0, 
Wye ey, aca = 164.0) 


Smith Predictor Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Raver CRUNMEMSL P23 se KG = 4 2.17.0, Ty = Ao/ 2.0, 


Ugg = 277.00, Kic 2 0/47, oes SUIS S 


PI Algorithm Control of Bottom Composition 
for a +25% Step Change in Feed Flow Rate 
(Run E-PI16; KRte= = 0-05.5118;," Ky = 201 81.22.) 


PID Algorithm Control of Bottom Composition 
for a +25% Step Change in Feed Flow Rate 
(Run E-PID23; Ke = -0.503, Ky; = -0.125, 

Kp = -0.075) 


Dahlin Algorithm Control of Bottom Compo- 
sition for +25% Step Change in Feed Flow 
Rate: e(Ruiny E-DAH2 6i30 Ky, =) gi 4); Tp = 918.0, 


Tete 1870 £ a0 1123: 10.) 

Smith Predictor Control of Bottom Compo- 
sition for +25% Step Change in Feed Flow 
Reaitiel a(R ulna vE ~S1Pi2 23st eh, (=) 422407 0), Ep) S720 , 
Tg t= isv20 Kp BO 71 dey. “Kp 8-011 85) 
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PI Algorithm Control of Bottom Composition 
for a -1% Step Change in Set Point 


(Run E-P129; NeU=5-0,607.45ka9e 70.0959 


PrDeALoOnichmecontroi] ot Bottom Composition 
for a -1% Step Change in Set Point 

(Run E-PID27; Ke eer OCS este ki n=Feod 1174 

Kp = +0200) 


Dahlin Algorithm Control of Bottom Compo- 
sition for °-1% Step’ Change in Set Point 
(Run E-DAH28; K, = -3.14, Tp = 918.0, 

eh = 91878 08 -xR=RY 6450) 


Smith Predictor Control of Bottom Compo- 
Str on wore} 1 Zstepy Change? inoSetm Point 
PRUNME—SEP3CESK Seee enangey = O184 0% 
ee eye Ose k=) 0669 = Kie=, -0¢ 2099 


PI Algorithm Control of Bottom Composition 
for a +1% Step Change in Set Point 


(Run E-PI31; Ke = 0.607, Ky = -0.095) 


PID Algorithm Control of Bottom Composition 
for a +1% Step Change in Set Point 

VRUNDEY PIDZOS EIR, SOM Ohh? 397 KT aL 20.1198 

Ky = - 0.200) 


Dahlin Algorithm Control of Bottom Compo- 
Sit OtmenOne tl eeotep Change in Set Point 
(Run E-DAH30; K, = -3.14, Tp = 918.0, 

THE EON 8700 .-452 Deano) 
Smith Predictor Control of Bottom Compo- 
Sit One OGetheeotep Change in Set’ Point 
(Run E-SLP30; Ky = -3.14, Tp = 918.0, 

Ty = 187.0, Ke = -0.869, Ky = -0.209) 


Dahlin Algorithm Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Rate (Run E-DAHO1; K, = -5.63, Tp =01760-0, 
Tit Lon ghar ,- 35% ose by 


Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 


Rate (Run E-DAHO2; K. = -5.63, Tp =p7'50'. 05 
Tor 8 1979) Shee Goseby 


Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Rate (Run E-DAHO4; K, = -5.63, Tp = 750.0, 
ee 8 L872) 2 Ver 234.5) 
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Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Rate (Run E-DAHO5; Kp = -5.63, Tp = 750.0, 
Ty = 187.0, A = 120.0) 


Dahlin Algorithm Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Rate (Run E-DAHO8; K, = -5.06, Tp = 750.0, 
Thue Fey lols be164.b) 


Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 


Rate (Run E-DAHO9; K, = -4.56, Tp = 813.0, 
eS 1972 O83 B= 16400) 


Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 
Rate (Run E-DAH10; Ky G38 8/4 Ip) =1853.0, 


ieee) 0s b= 64. 0) 


Dahlin Algorithm Control of Bottom Compo- 
Sition tor -25% Step Change in Feed Flow 
Ratemqnunge-DAHII eK, 72° 23,49» Tp = 896.0, 


p 
Th 219977 OM E"1 64105 


Dahlin Algorithm Control of Bottom Compo- 
sition for -25% Step Change in Feed Flow 


RatemoRUn E-DAHL2s. Kk t= 623.14), Tp = 918.0, 
UF) Gee AUS at age 164.0) 


Smith Predictor Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Ra CemURUMe= SE RODS Ks = 2314), Tp = 918.0, 
Tq = ie hy ceed Oh Kc = -0.830, Kt = -0.209) 


Smith Predictor Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 
Ra temuRUnebesLl PLO si Ke to 2301 4. leet S 30, 
Vieem So /Use Ku SO B69y Kl= 00.4209) 


Smith Predictor Control of Bottom Compo- 
Sition for -25% Step Change in Feed Flow 


Ra Ce URUME = SUP L9 eK core oe ae Tp = 918.0, 
Wee Te 1G = enn, fe 8. STS 


PID Algorithm Control of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 
(Run E-PIDO5; Ke = -0.636, Ky = -0.200, 

Ky = -0.666) 


PID Algorithm Control of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 
CRUMBe=b DOI Rhee =! 057015 5Ky) = 0.105; 
Kjos 0. 333) 
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PID Algorithm Control of Bottom Composition 


for a -25% Step Change in Feed 
(Run E-PID13; Kppar 05035) Ky 


Ky senor200) 


Flow Rate 
SO. LUOn 


PIDSALGoeithmeControl,of Bottom Compostbion 


for a -25% Step Change in Feed 
(Run E-PID15; Kee= 30.659 7K] 
Ky = -0.170) 


PIMA GonrTchnm,control of Bottom 
for a -25% Step Change in Feed 
(Run E-PI14; Ke = -0.551, K] = 


Puen ooGithm Control ofr Bottom 
fora =25% Step Change in Feed 
(Run E-PI18; Ke = -9.559, on 


PIVALgonithm Control of Bottom 
for a -25% Step Change in Feed 


(Run E-P120; Ke =)-0.551, Ky = 


PI Algorithm—Control of Bottom 
for a -25% Step Change in Feed 
VR UNME Cree 2geK ea e0'.578, KI = 


Flow Rate 
At) la 


Composition 
Flow Rate 
SOc.) 
Composition 
Flow Rate 
Ore t) 


Composition 
Flow Rate 
-0.106) 
Composition 
Flow Rate 
-0.100) 
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List of Symbols 


Symbols for Control Algorithms 





Symbol Description 


e: GOetinrcientsrinvetica numerator polynomial of the 
digitale concrouler transter function D(z) 

uF coeificients in the denominator polynomial of the 
Gigutaltcontroiler*transter function D(z) 

:; coefficients, unique to particular equation 

G controlled variable 

D GEanstens.unctionyofktaadigitaljcontrokler 

e error = set point - measurement reading 

E Caron iN sSeOnez domains 

£ coefficients, unique to particular equation 

G process transfer function 

Gy process transfer function 

Ge oh transfer function of a zero-order-hold device 

H transfer function of a measurement device 

k Sample index 

K transfer function of the close loop response 

Ke proportional *sconstdnt ®inesplenror ePlDAScontral 
algorithm 

Kp derivative constant in PID control algorithm 

KT integral constant in Pl or PED control algorithm 

Ky process gain 

m order of a polynomial 

M measurement time delay expressed as a number of 
Sample intervals, M = T/T 

n order of a polynomial 

N total system time delays in number of sample 
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S 


At 


process time delay in number of sample intervals, 
Reet Thy 


set point in s or z domains 
Laplace domain variable 
time domain variable 


Sample time 


closed loop total system time-delay 
process time delay 
derivative time constant in PID control algorithm 


gral time constant in PI or PID control algo- 
m 


measurement time deiay 


process time constant 


contr olmtere tout put 
CONntCrOller sOuLDUL <In Sor z domain 


z-domain variable 


Description 


System matrix, defined by component balance equa- 
tions 


System matrix, defined by mass and energy balance 
equations 


integration time interval (s) 
Murphree vapour efficiency 
feed mass flow rate (g/s) 
liquid enthalpy (Kd/g) 

vapour enthalpy (Kd/g) 


liquid mass flow rate (g/s) 
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Q heat input (Kd/s) 

QL heat loss (Kd/s) 

Qp reboiler duty (Kd/s) 

R reflux flow rate (g/s) 

S Side stream mass flow rate (g/s) 

V vapour mass flow rate (g/s) 

WT ' liquid mass holdup in a stage (g) 
x liquid composition (wt.% methanol ) 
y vapour composition (wt.% methanol ) 


Subscripts 


7] initial value 
i? final value 
n stage number 


Superscripts 


k variable value at the kt" integration cycle 


% equilibrium composition 
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1. Introduction 


1.1 Objectives 

The availability of low cost digital computers with a 
high degree of reliability, particularly the new generation 
of process computers that are now available, is having a 
Significant impact on computer utilization in the process 
industries. For control engineers, it means that more 
advanced control schemes and control algorithms can _ be 
implemented by the process computers. This project deals 
with an evaluation of different digital control algorithms, 
namely the deadbeat algorithm, the Dahlin algorithm, the 
discrete Smith predictor, and comparison with the discrete 
proportional-integral (PI) and proportional-integral- 
derivative (PID) algorithms. 

In the implementation of digital control algorithms, 
some of the primary concerns in selecting a general digital 
control algorithm are: 

1) the algorithm should easily be a standard soft- 
Warentbhockesof elaediinest, digital. control (DDG) 
package; 

2) the algorithm should be easily programmed; 

3) the computational efficiency of the algorithm 
should be high; 

4) the output of the algorithm should be calculated 
only based on the present error, past errors and 


past outputs; and 
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5) the algorithm should be general enough that 
different control algorithms can easily adapt to 
Ls ihorm. 
In evaluation of the performance of different control 
algorithms, the primary concerns are: 
1) the performance under regulatory control; 
2) the ability of the algorithm to compensate for 
process time delay; 
Smecne ssimplaicltyeot yon-linentuning: and 
4) the influence of model errors and nonlinearity on 
the performance of the control algorithm. 

The experimental work of this project was carried out 
Oumeecne pi lotvemscane sbinary “distillation column in. the 
Department of Chemical Engineering. Control of the distilla- 
mrOnee COLUMN Iss) ditticult. .sinces it exhibits ‘nonlinear 
behaviour as well as time delay. These process characteris- 
tics plus the increased time delay due to the use of a gas 
chromatograph for bottom composition analysis make _ the 
column a realistic system for the evaluation of the perform- 


ance of different control algorithms. 


1.2 Thesis Organization 

This thesis iS organized in the following manner: In 
Ghapters 2. teasabrietmelitecaturesisurveyeson digital control 
algorithms is presented; the development of the various 


digital control algorithms is given in Chapter 3; Chapter 4 


describes the development of a nonlinear binary distillation 
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column model; the results of the digital simulation studies 


On the performance of different control algorithms are 
presented in Chapter 5; implementation of the digital con- 
trol algorithms is described in Chapter 6, as well as the 


Results eoferthe erexpenmimental rewaluatiion:s ofsothen control 


algorithms; finally, the overall conclusions and recommen- 


dations from this work are summarized in Chapter 7. 
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2. Literature Survey on Digital Control Algorithms 


2.1 Introduction 

With the flexibility provided by digital process com- 
puters, control engineers can now implement not only more 
sophisticated control strategies, such as different levels 
of supervisory control for plant-wide optimization, but use 
advanced control algorithms for low level control to achieve 
better process performance. In this chapter, a brief litera- 
eunesssurvey of*'tne ‘digital’ 4Ycontrol’> algorithms}? for’ ‘single 
input single output feedback control systems, which are 
readily available for industrial applications will be dis- 
cussed. The format of the discussion is based on the evolu- 
tion of digital control applications in the process indus- 
tries. The survey covers developments from the use of the 
digital equivalent of conventional control algorithms to the 
present involvement of designing digital control algorithms 
using sampled-data control theory. Tuning techniques for the 
algorithms cited in this chapter will also be discussed with 
emphasis on the tuning procedures for regulatory control 


applications. 


2.2 Digital PI and PID Algorithms 








The digital proportional-integral (PI) and proportional 
-integral-derivative (PID) control algorithms were the ear- 


liest digital control algorithms applied to control applica- 


tions in process industries when the digital process comput- 
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eansyiweregermtnuaciced font spl antaycontmolt.. Duringtethe Bperiod 
from the early 60's to the mid 70's, methods of obtaining 
the digital equivalent of the conventional analog PI and PID 
controllers and the tuning techniques for these two digital 


control algorithms received the most attention. 


2.2.1 Positional and Incremental Forms 

There are two basic forms, namely the positional form 
and the incremental form (also known as the velocity form), 
of digital PI and PID control algorithms. The most commonly 


used positional form of an ideal digital PID control algo- 


rithm is described by the following equation: 


li TK Tp 
z Cat 
i=1 aT. 
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Mi 
Its counterpart, the incremental form is obtained by 


taking the difference between Equation 2.2.1 and the equa- 


tion for the previous sample (k-1). The resulting expression 


is 
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In general, the incremental form of the ideal digital 
PID controller is preferred rather than the positional form 
because of its "“bumpless" transfer and anti-reset windup 
features [Bibbero (1977), Bristol (1977), Desphande and Ash 
€2981.) weSmirthi (dl 972) Je 


Different approximation techniques for the _ integral 
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and/or derivative term(s) have been suggested [Bibbero 
C197 7s Wirt neome ye, Verbruggen, etsaitw(1975)) ctcieimprove 
the performance of the algorithm. A summary of the different 
forms of digital PID control algorithms, in the z-domain, 
that result from using different discrete integration tech- 
niques is presented by Verbruggen et al. (1975). Although it 
has been claimed that using different approximation tech- 
niques, such as the trapezoidal rule for integration [Bibbe- 
ro (1977), Smith (1972)], the four point central difference 
formula for the derivative [Bibbero (1977)]} would improve 
the operations of the algorithm, no comparison of these 
"improved" algorithms with the ideal digital PID (cf. Equa- 
tions 2.2.1 and 2.2.2) has been undertaken to the author's 
Knowledge. 

Another variation from the ideal PID algorithm is the 
Change of the feedback variable to the derivative term. 
Smith (1972) commented that "there is no real advantage in 
allowing the derivative mode to nan ae a change in set point 
Signal", therefore, the feedback variable for the derivative 
term, as suggested by Smith, should be the measurement 
sequence rather than the error sequence. The advantage of 
the resultant algorithm is the elimination of a "kick" from 
the derivative action due to a set point change and the 
algorithm is essentially the same as the original form when 
it is applied to regulatory applications. Other considera- 
tions regarding implementation of digital PI or PID algo- 


rithms have been discussed in detail by Bristol (1977) and 
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Smith (1972), but some of the material, particularly that 
dealing with hardware limitations, is now out of date. 
Despite the existence of various forms of digital PI 
and PID algorithms with practical enhancements, the digital 
PI and PID algorithms still most frequently used for compar- 
ison with the performance of control behaviour achieved 
using other control algorithms are those given by Equations 


Saath anid (24.123.3% 


eee. 2flunings Digital+Pr<and) PID) Algorithms 

The main advantage of using digital PI and PID control- 
lers over the other digital control algorithms is that they 
are general purpose single loop controllers. However, tine 
wuning ther Piytand, RD controllenssto-yield optimal) control 
performance is not a simple task because of the necessary 
trial and error procedure. With the appearance of low cost 
digital computers, finding the optimal controller settings 
uSing a systematic approach can be performed through digital 
Simulations. Lopez et al. (1967) proposed that’ error- 
integrals of the closed loop response of a control loop 
Should be used as the tuning criteria for analog PI and PID 
controllers rather than the quarter decay ratio criterion as 
proposed by Ziegler and Nichols (1942). Three types. of 
error-integral, namely the integral of squared error (ISE), 
the integral of the absolute value of the error (IAE) and 
the integral of time multiplied by the absolute value of the 


error (ITAE), were suggested as the criteria for adjusting 
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the controller settings of the analog PI and PID control- 
bers. The= characteristics’ of the’ different’ integrals and 
their applications were also discussed by Corripio et al. 
CP97s)eincaliiater publ icationeeithnough Simulation «studies; 
Lopez et al. (1967) were able to relate the parameters of a 
first order plus time delay transfer function model with the 
controller settings for analog PI and PID controllers which 
could be used to minimize the error-inteygrals and optimize 
the control performance. The correlation between the trans- 
fer function parameters and the controller constants were 
Givens int equativonetorm.PForidigitaloplrandePpiDécontrobilers,; 
Lopez et al. (1969), using the same approach for the case of 
the analog PI and PID controllers, generated a set of tuning 
graphs correlating the model parameters and controller set- 
tings which could minimize the error-integrals. A similar 
approach has also been used by Roberts (1976). Unfortunate- 
ly, convenient equations were not reduced from these graphs 
due to the additional parameter, the sampling time, present 
in the control algorithms. For a control loop in which the 
Sampling tet ect) tisstinsignifticdnt panies thelratiomof thepsamns 
pling time to the most dominant process time constant is 
Smaller than the ratio of the process time delay to the most 
dominant process time constant, correlations for the analog 
PI and PID controllers may be extended to the digital PI and 
PID controllers. Smith (1972) suggested the following proce- 
dure: add one half of the sampling time to the process time 


delay to compensate for the sampling effect and then apply 
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the correlations for the analog PI and PID controllers to 
obtain the controller settings for the digital PI and PID 
controllers. Other tuning approaches that do not use error- 
integral criteria have been employed. Chiu et al. (1973c) 
Suggested the controller settings for PI and PID controllers 
be calculated using Dahlin's analogy that the desired closed 
loop response be specified and the parameters in the control 
algorithms be calculated from the closed loop characteristic 
equation (the Dahlin algorithm will be discussed in a later 
Section). Vergruggen et al. (1975) summarized the coeffi- 
Evenes (Ofme hetidiscreteeP lrandsepIDicontrollers as ae function 
of the ultimate gain and ultimate period of oscillation 
which are the Ziegler-Nichols closed loop criteria (1942), 
with adjustments for the sampled-data control system. Coef- 
fercaents “or thet dascretes Pl sande PID! tcontrolic algorithms 
according to the Haaiman adjustments for analog PI and PID 


controllers have been reported by Verbruggen et al. (1975). 


2.2.3 Choice of Sampling Time 

In the digital PI and PID control algorithm, sampling 
time is an important factor which can affect the performance 
of the control algorithms significantly. Choosing the appro- 
priate combination of the sampling time and the controller 
settings is a difficult task. In industrial practice, the 
sampling time for a control loop is chosen primarily accord- 
ing to the process conditions, e.g. one second sampling time 


for a flow loop [Smith (1972)], as opposed to the rigorous 
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analysis which is based on the stability of the control 
loop. The choice of sampling time from a practical point of 
view has been discussed by Smith (1972). Verbruggen et al. 
(1975) present a simple correlation for relating the choice 


of sampling time to the process time constant. 


2.3 Discrete Smith Predictor 

Time delays, such as transportation delay, measurement 
delay, are inherent in many chemical processes. When the 
Cime delay usisignificanthtim a: control: loop, ise.sctheeratio 
of the time delay to the process time constant is relatively 
large, the control pertormance that is achieved using a PI 
Or PID controller deteriorates. One of the special control 
schemes for controlling time delay processes that’ has 
received the most attention is the Smith predictor [Smith 
(1957,1959)]. The basic philosophy behind the development of 
the Smith predictor is to introduce an element into the con- 
trol loop in such a way that the time delay element does not 
exist in the closed loop characteristic equation. Unfortun- 
ately because of the hardware implementation problem this 
scheme was not feasible until computers came into use in 
process control applications. 

Experimental evaluation [Alevisakis and Seborg (1974), 
BUCK ve Vutl 900 ee DOSS (19/74) m0osss andy Moore. (1973), “EBupfer 
and Oglesby (1962, 1961), Meyer et al. (1978), Prasad and 


Krishnaswamy (1975)] and simulation studies [Meyer et al. 


(1976), Nielsen (1969), Schleck and Hanesian (1978), Smith 
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and Groves (1973)] have demonstrated that for controlling 
time delay processes, the Smith Predictor can improve the 
control performance significantly as opposed to that achiev- 
ed by using the conventional feedback controllers. Sensiti- 
vity analysis of modelling errors in the Smith predictor has 
also received considerable attention [Buckley (1960), Eisen- 
berg (1967), Garland and Marshall (1975,1974), Marshall 
(1974)]. For practical implementation, Meyer et al. (1978) 
presented the derivation of the Smith predictor using time 
domain analysis while Deshpande and Ash (1981) developed the 
algorithm in the z-domain. Both approaches are well docu- 
mented and the final equations are essentially equivalent. 
The discrete Smith predictor, in both cases, is in a form 
ready for programming, however, it should be noted that nei- 
ther of these two approaches reduces the Smith predictor 
into the form of a single equation. Consequently, some ad 
hoc procedures are required for implementation of these 
equations with any general digital control package. Thesen 
(1980) presented an alternate derivation of the discrete 
Smith predictor, however, he concluded that the alternate 
form, which is more difficult to derive, does not yield any 
better control performance than that obtained using the dis- 
crete Smith predictor previously developed by Meyer et al. 
(1978). 

For tuning the Smith predictor, the common approach is 
to assume that the model used in the algorithm is accurate 


So exact cancellation of the time delay term occurs. Conse- 
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quently the control constants in the PI controllers can be 
designed based on the effective closed loop characteristic 
equation or any correlations for PI control of processes not 
considered to contain time delays. In the experimental eva- 
luation of the Smith predictor, Meyer et al. (1978) used the 
correlations for PI controllers suggested by Moore (1969) to 
Obtain the initial tuning constants and the controller was 


then fine tuned for optimal performance. 


2.4 Sampled-Data Control Algorithms 
2.4.1 Direct Synthesis Method 


While the digital equivalent of conventional analog 
controllers were implemented quite successfully on digital 
process control computers, digital control algorithms de- 
Signed from an entirely different approach was under inves- 
tigation. Koppel (1966) claimed that "there is no real rea- 
son to retain the proportional-integral or _ proportional- 
integral-derivative forms in direct digital compensation. 
These forms were developed for analog instruments. The 
greater flexibility made possible by use of compensation 
Should be utilized in designing digital algorithms". The 
particular design method for digital or sampled-data control 
systems suggested by Koppel was the direct synthesis method 
developed by Ragazzini and Franklin (1958). 

In «the design»s:procedure, using the direct synthesis 
method, ‘requires that the transf_r function of the process 


be known, then the controller can be obtained through the 
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manipulation of the overall transfer function providing the 
transfer function of the desired closed loop behavior is 
Specified. Applications of this design procedure have been 
peesentedsnnsivarying degrees, offedetail? byscBadavase()981 ); 
SuiuUSClOPN) BEChiO ety aleonngo79c) seCondon! andtSmith (C1977b))} 
Dahlin (1968a), Deshpande and Ash (1981), Koppel (1966), and 
Lane (1970). The major advantages of using this design tech- 
muque sarenarcsansimplicuuy ‘and? 7Tlexibi lity yttive.eddififerent 
control algorithms can be synthesized depending on the spe- 
eurications of® thea transfer | functionarfor? ithececlosed sloop 
response. 

It should be also noted that since the control algo- 
rithms designed using the direct synthesis method begins 
with a discrete model established for a particular sample 
mcenval, the -performance of the control algorithms is 
affected by whether the discrete transfer function model, 
which is sample interval dependent, represents the actual 
DeoGess tafaimtyesiLvuny 161956); ekuo) eGb9/74y, CRagazzinig§gand 
Franklin (1958) ] 


2.4.2 Deadbeat Control Algorithm 

The deadbeat control algorithm is synthesized using the 
direct synthesis method by specifying the transfer function 
for the closed loop response to a step change in set point 
as K(z) = Zen ale where n is the number of sample intervals 


due to the time delay of the process and 1 designates the 


time delay due to sampling. The physical interpretation of 
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thew abovertranstereTrunction Wse*that’saesthe n+] th Sampling 
bastantrsatternethe bintroduction, ofi@atrstep echange*iini set 
point, the deadbeat control algorithm will bring the con- 
trolled variable to the desired set point and will remain at 
the set point for all subsequent sampling instants. Although 
Simulation results [Chiu (1971), #4Koppel (196 6. eo taSm neh 
(1972)] were promising, however experimental results report- 
ed by Uronen and Yliniemi (1977) have demonstrated that the 
deadbeat control algorithm did not provide as satisfactory 
performance as did the digital PI and PID control alygo- 
rithms. Application of the deadbeat control algorithm in the 
process “industries is*also not ‘popular primarily due to the 
fact that it lacks convenient on-line tuning parameter(s) to 


adjust for changes in process conditions. 


2.4.3 Dahlin Control Algorithm 

A digital control algorithm which is gradually gaining 
industrial acceptance is that developed by Dahlin (1968a). 
This single input single output algorithm was also proposed 
for the control of multivariable systems when the control 
system is completely decoupled [Dahlin (1968b)], however 
this multivariable control scheme did not receive much at- 
tention probably due to the fact that the frequency domain 
design technique for multivariable control systems was dom- 
inant at the period. In the subsequent discussion, only the 
Single input single output version of the Dahlin algorithm 


will be presented. 
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The Dahlin algorithm is also synthesized via the direct 
Synthesasa method. alnmethis) casejythe sransfern.function. for 
the closed loop response to a step change in set point is 
Specified to be first order plus time delay. With the mani- 
pulation of the time constant, A in the closed loop response 
reaMsnerwabunct ion se .desired.vpentormance..can,. be, obtained. 
Also, manipulation of this variable provides on-line tuning 
Ofercherscontrol calgonttnm -for ditferent. process, conditions. 
The ringing free form.of the Dahlin algorithm can be obtain- 
ed by removing poles that are close to the unit circle [Chiu 
Pie ailae@h973a)5 Dahlin) s(1968a) |.» It.eshould,.bes noted that 
this ringing free form of the algorithm has the same struc- 
eer, cthesdigital Epil ons Pidi.control algorithmy.if gagfirst 
Or second order plus time delay process transfer function is 
used respectively. 

SamuUlidiioOnme wesutusiBkadavas (1981), Chiu.(19/1), Chiu 
et al (1973a), Condon and Siiishelie AMO 7 1/0), ee t/a) Dele em oman 
(1972)] have shown that the performance of the Dahlin algo- 
rithm for both servo and regulatory control operations is 
comparable to that achieved by using the digital PID control 
algorithm. For a control system with a significant time de- 
Layee the iper fonmancel cusingabhesDahlin algorithm. is. round to 
be superior to the performance that is possible using a di- 
gital PID control algorithm. On the other hand, experimental 
results obtained by Uronen and Yliniemi (1977), showed that 
the performance of the Dahlin algorithm is rather sluggish 


when it is used for regulatory control. These conclusions 
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however are subject to the specific control systems being 
considered since a general analysis on the performance of 
the control algorithm, from a theorectical approach, to the 
author's knowledge does not exist. Sensitivity analysis to 
demonstrate the effect of inaccurate model parameters used 
in the Dahlin algorithm have been undertaken by Badavas 
(1981) and Condon and Smith (1977b) using digital simula- 
tion. From the simulated results, Condon and Smith observed 
that the performance of the Dahlin algorithm for servo oper- 
ation is improved noticeably if the process time constant or 
the process time delay is over- estimated. However, an inap- 
propriate value of the process time delay will cause the 
control system to be unstable. The sensitivity analysis 
study by Badavas reported stability limits which correlate 
with the parameters in the process transfer tunction from 
frequency domain considerations. 

On-line tuning the Dahlin algorithm is quite straight- 
forward. After obtaining the best estimated process model 
parameters, the initial value of A} is set to be the same in 
value as the process time constant. The desired performance 
can then be achieved by adjusting A according to observa- 


tions of the closed loop response. 


2.4.4 Other Sampled-Data Control Algorithms 





An algorithm which can be considered as a variation of 


the Dahlin algorithm is that reported by Higham (1968) and 


Jones (1978). The primary difference between the Dahlin 
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algorithm and the one of Higham is that the incremental out- 
put is calculated in Higham's algorithm as opposed to posi- 
Civonavesoutput’ ebeing® calculatedttinm thet Dahlin algorithm. 
Tsing et al. (1979) proposed another variation of the Dahlin 
algorithm based on Higham's analogy; instead of uSing a 
transfer function to represent the process, Tsing and co- 
workers suggested the use of the actual data of the process 
response to a set point change in their algorithm. The 
Pesultant= algorithm ts’ rr*terms’*of the’ present’ ‘error term 
and the incremental output sequence. Both the proposed alyo- 
rithms provide a single on-line tuning parameter similar to 
that in the Dahlin algorithm. 

The Kalman algorithm [Chiu (1971), Chiu et al. (1973b), 
Smith (1972), Uronen and Yliniemi (1977)] is designed on the 
basis of a specification of both the responses of the mani- 
pulated and controlled variables. Simulation studies by Chiu 
(1971), Chiu et al (1973b), Smith (1972) and experimental 
evaluation by Uronene and Yliniemi (1977) have shown that 
the Kalman algorithm will also provide satisfactory perform- 
ance. However, the Kalman aigorithm suffers form the same 
disadvantage as the deadbeat algorithm, namely that it is an 
algorithm that does not contain parameters suitable for 
On-line tuning. 

In the design procedure for the different control algo- 
rithms that haven been discussed, a load disturbance to the 
System is not considered although some algorithms do have an 


on-line tuning parameter to adjust for changes in process 
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conditions. Mutharasan et al. (1978) proposed an algorithm 
which can compensate for both set point and load changes. 
The design procedure involves determination of values for 
the manipulated variable to achieve a desired response in a 
manner which is analogous to the direct synthesis method. 
The control alyorithm is also a simple recursive equation of 
input and output sequence. Although the algorithm is a 
Steady state feedforward type of control algorithm, it 
Should be possible to easily extend to a predictive type 


control scheme. 


2.5 Future Trends 

From the preceding discussion, it can be concluded that 
a control algorithm designed with at least some knowledge of 
the process, ie. a mathematical model which represents the 
dynamics of the process to be controlled is included in the 
design procedure, can provide better performance than the 
conventional Bea por] PD acontro|l fal gorithms.¢ On” thedsother 
hand, it is also understood that the use of an inaccurate 
model in the control algorithm will cause the performance to 
deteriorate significantly as well as tend to cause system 
instability. For these reasons, adaptive control schemes 
become very attractive to control engineers. 

One of the adaptive control algorithms which has 
received considerable attention in recent years is the 
self-tuning controller. Application of self-tuning control- 


ler on distillation column can be found in Lieusion (1980). 
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3. Digital Control Algorithms 


oe) -[ntroguckion 

In recent years, due to significant hardware improve- 
ments, there has been a rapid growth in the use of digital 
computers for process control. Although control engineers 
Can.now implement without difficulty more advanced control 
senemes oretcontrolotal gorithms; imostftota the control loops sin 
industry still employ only conventional single variable 
meedbackycontroleusing Pil “or, PdDscontrol.-algoniithms., Part -of 
the reason for this situation arises from the fact that many 
of the advanced control algorithms designed for continuous 
control systems often require an ad hoc approach for a digi- 
tal system implementation. 

In this chapter, a general purpose digital control 


algorithm of the form 





U(z) ao re yy rages Doers ame rele + an zeae 
at , a + + b -m+l 4 p 77m 
EZ) bg + byz so m-12 m4 
(eSeue cua) 


will be used for presentation of the different algorithms. 
Advantages of using the above algorithm for computer 
process control are: 
1) it can be easily programmed; 
2) the computational efficiency is high, because it 


involves only simple arithmetic operations; 


SecnecercCONurolmoucpul Ulatsmecalculated based’ on ‘the 


present error, past errors and past outputs; and 
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4) almost all control laws can be expressed using 
Chins stOrneor algorntuthm. 

Since a system operated under computer control is 
sampled-data in nature, the mathematical derivation involves 
manipulation in the z-domain. The direct synthesis method is 
chosen in this work for designing digital control algorithms 
because of its simplicity. This technique will be presented 
in the next section. The digital form of the deadbeat and 
Dahlin algorithms, which are derived by this technique are 
discussed in Sections 3.3 and 3.4 respectively. Since the 
Smith predictor has gained acceptance in industry for the 
control of processes containing time delays, the discrete 
form of this well known algorithm is given in Section 3.5. 
Peune laste Ssecuton Of Unis Chapter, Section 3.6, the incre- 


mental form of the PID algorithm is presented. 


3.2 Direct Synthesis Method 


—_. 





The direct synthesis method is a simple and straight- 
forward design procedure. This method can be applied to 
either continuous or sampled-data control systems. For the 
control system shown in) Figure 3.2.1) the closed lop) trans- 
fer function in either the s-domain or the z-domain is des- 


cribed by the following equation: 


G DG 
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Assuming that the transfer functions G@ and K(=C/R) are 
known, then the expression for the controller D is given by 


the following equation: 


U 1 
eee (2b) 
E Gi Valea 


Pemicaneoe Seen from Equatijon 3.2 2 that different con- 
trol algorithms can be established depending upon the choice 
of K as selected by the designer for a particular applica- 
tion. 

The direct synthesis method can also be extended to 
control systems which do not have unity feedback, for in- 
Suanece,atne contra! system shown inyFrgure 3.2.2 contains a 
measurement device which has a transfer function other than 
unity. 

PhemeclasedssOOpemethansier | function ofS this system 


expressed in z-domain representation is: 





a ae D(z) G(z) (ness 

R(z) ee DiiZ1)) Giz es Haz, ) 

where G(z) = GAB onee:° Gp(s)] (cea) 
H(z) = GLH(s)] (3.2.5) 


Sor it; follows athnat the. transfer function “oT the digital 


controller D(z) derived by the direct synthesis method is: 











U(z) 1 K Cz) 
D(z) = = ——$—$—————— (or. 6) 
Eetez,) G(z) TetaeKe(e7 ali Z2) 
Giz) 
where («2 i= (Serres) 
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3.3 Deadbeat Control Algorithm 
Deadbeat or minimal response is one which satisfies the 
following criteria (28): 
1) The system must have zero steady-state error at 
Chee SanDLLOGeauiS uous on the specified tnput 
Cesc ssilignads, 
2) The rise time should equal a minimum number of 
Sampling periods. 
3) The settling time, measured at the sampling 
instants, should be finite. 
Cieza Zeeandi Gz) / RUZ) emust. alii bes physical ly 
realizable. 
FOr a minimal phase system, the deadbeat response to a 


Seeominput-in Set point is (28): 





K(z) = = Sit (OF seit) 
R(z) 


However for a system with a process time delay which is 
greater than or equal to P sampling intervals but less than 


P+1 sampling intervals, the response is: 


ers ee (orga) 
R(z) 


For the sampled-data control system illustrated in 
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Goh(s) = ——— (Beach) 


H(s) = exp(-T 5) oro 10) 


where Ky = process gain 


process time constant 

Ty = process time delay 

m = measurement time delay 

T = sampling time 
werescrcucton = OT sEquations 3.3.3 "“and=3.3.4° Into’ Equation 
See.d, yields 


ela) @ Sot (353.65) 


oe eae at 
fy = exp(-T/Tp) 
2 = exp(-mpT/Tp) 
M5 Sele ae ust, P 
P = ahi (integer) 
For a measurement device with a transfer function given by 
Equal ion 3.3.5, LnGmz-domaingrepresentation 1S: 
HZ )e =m 
where M = Ti/T (integer) 
Bor, othe sitransrer = etunctions  detinedm Dy™ Equat1onsy 3.2.6, 
OeisticuOeoeOmandes 25. suseH Of stoner hzects synthests ‘method 


gives the deadbeat control algorithm listed in Table 3.3.1. 
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Ika DilWe met aes 1h 
DEADBEAT CONTROL ALGORITHM 


U(z) dg + ayia 
= secede CEREMERTS EET ie Na (seo) 
E(z) Dy + byzve BHD PL Biot Dayo zOM 
1 = 
= : Lage(k) + aye(k-1) - byu(k-1) - byyyu(k-N-1) 
0 
- Byegu(k-N 2)] (Soa 
where ay = ] 
pr eecrec 01 
Do =a CH 
DF — Kpco 
Py+r = Kped 


1 = exp(-T/Tp) 

2 = exp(-mpT/Tp) 

1 An ee itl 

Fede Padua din Ml 

D = eid/ sy P 
= lial sCinteger) 
= t/a (integer) 
= P +M 


as 
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3.4 Dahlin Algorithm 





For most industrial applications, deadbeat response is 
ideal but usually not realistic because of process condi- 
tions. Furthermore, the deadbeat algorithm suffers from the 
fact that it does not have any provision for convenient on- 
line tuning of parameters. An algorithm that contains a sin- 
gle parameter that can be conveniently adjusted on-line is 
that of Dahlin (15). Dahlin proposed that the closed loop 
mesponse of a control system for a step change in set point 
Should behave like a continuous first order plus time delay 
mEOeeSS With the time constant of the closed loop transfer 
function being the on-line tuning parameter. The closed loop 


system response, for a System controlled by the Dahlin algo- 


rithms, is given Dy 





C(s) exp(-T.s 
Armia?! (gaara 
R(s) NS et oe! 
in the s-domain, and by 
exp(-T S 1 
fo any 
Caz) Xd sAtL 1 S 
ee ee (Sa4i323) 


1 
Ce 
S 
in the z-domain. Bin sEquatironm 34-2, Ve is the time delay in 
the closed loop response and X is the time constant of the 
closed loop response which is the tuning parameter. 
Since the controller must be realizable and must not 
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Equation 3.4.2 must be chosen no smaller than the value Ty 
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Pee CUatOnms scr Omertnes Limitinds case, Equation 3.4.2 


becomes: 
K ie(iey “ec Bae) zrbol 
Ri ee ee eer ean Ce ore a pgeaesy 
(1 - pee) 
where Cc, =" &3 f5 
a= to Rah 
fy = exp(-T/a) 
fr, = *expiiem |r/2) 
p 
= 1- gif | rer 
p 
P = Td/T 


FOr the sampled-data controlc system considered in Section 
S-2, a¢pplyinguthe dimect ssynthesis method using the transfer 
meno Ols eecet ned B Dy Pebquations. 3.2.77,  3.3-6, .3.3./7- and 
wees seune “wahlim: e€0ntro! algorithm can be obtained. The 
Dahlin algorithm derived by this procedure is given in 


fabye  3.4.1- 


Laie 3 4" 1 


DAHLIN ALGORITHM 
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JL SV te rae Rojeckw rate k-1)) taoe(k-2)) = byu(k-1) = 
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i) Sy Se aera 


Kpc3 

Kp(cq - ¢3¢)) 
-Kpc4fy 

-KpC3c1 

-Ky(Cqcy tac coy) 


Kpc4c2 

Pus aay 

Wa nel 

l= 4 

fq - f3 
exp(-T/i) 
exp(-m Tia) 


p 
exp(-T/Tp) 
exp(-m /Tp) 
1 - t/t + P 
Tq/T 
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See iscratvesomthaPpredictor 





The Smith predicton is. well. known: for its ability to 
improve the control performance of processes that exhibit a 
Significant time delay. The derivation of the Smith predic- 
tor has. been given by numerous different workers (cf Chapter 
Sieeoue Only fom continuous systems. ‘Inithis section, the 
discrete form of the Smith predictor is derived. The block 
diagram representation of the Smith predictor for continuous 
enae sampled-data control systems is given in Figures 3.5.1 
anon. .5.2 respectively. 

in order to simplify the derivation of }the ‘discrete 
eentne predictors it is convenient! to use a different block 
ivagran than that given in Firqure™seo.2. By block diagram 
reaquetion,: the block diagram shown in) Figure;3.5.2 can. be 
mdgar ied to tna shown in Figure 3.5.3. 

Prom ihidune 32523, it follows that jthe expression for 
the discrete Smith predictor is 
D(z) = 2. 2 /Ne 2 ee (G0) 

ene z) 20> (-z) 
with az) the discrete PID algorithm, given by 
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aque ajyz7- + ae 


"Ge eae aan CS sia) 
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ay = -(Kr + 2Kp) 

a> = Kp 
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and Do(z), the time delay compensation part of the Smith 


predictor, expressed as 


GiGp(s) (1 - expl-(Tg + Tm)s] —)] 
Do(z) = Ta iki, ea ee LS (orb. 3) 
Glinae 
where Go(s) is [Kp/(Tds + yl), athe process transfer function 
(cf Equation 3.3.3) without the time delay element. Substi- 
Pim LoOneOt- CihlSe expnessi on ,intorsEquation i3.5.3°,ablows the 


time delay compensation part of the Smith predictor to be 


written as 


egieeeer “\ rere es (Cece ze!) zach] 
p 1 1 2 


where f, = exp( T/Tp) 
9 = exp(-MpT/Tp) 
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bReozeto y Degas + Dyygqz 
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(Shoo 5) 


= 26 [age(k) + aye(k-1) + age(k-2) + a3e(k-3) - 


bju(k-1) - bou(k-2) - b3u(k-3) - by, yu(k-N-1) - 


Dysou(k-N 2) - byygu(k-N 3) = byyqu(k-N 4)] 


where ay 


a0 

a, - fiag 
a2 - fia} 
-fia9 


Caen a eectan  [) 

-(K] + 2Kp) 

Kp 

1 

Spel ee A er) 
ncn: Soares aT 
etal ET) 
~Kpagcl 

-Kp(a,cq + agcCo) 
-Kp(a2c, + ajc?) 
Kpa2c2 
exp(-T/Tp) 
pute ep) 


(325.0) 
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f.0) Discrete PID Algorithm 





Proportional-integral (PI) and proportional-integral- 
derivative (PID) control algorithms are the most commonly 
used control algorithms in industry, however, there are 
numerous digital equivalent forms of these algorithms. The 
different forms result basically from different numerical 
approximations of the integral and derivative terms. Never- 
theless, the discrete algorithms can be divided into two 
main categories, namely the positional form and the incre- 
mental form (also known as the velocity form). In this work, 
the incremental form of the ideal PID algorithm is chosen 
Porm. the trol lowind reasons: 

byte do=senot require a Dias term as mn) ‘the sposae 
GrOnals form: 

2) it allows for “bumpless" transfer from manual to 
aucOlba GiCe CON GEOne 

3) it provides for anti-reset windup; and 
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The continuous ideal PID algorithm can be written as 


1 de 


u(t) = Kefe(t) + ofie( tye dt eT p snd een Oak) 





I 

where U.. = steady state output or bias 
The digital equivalent of the above equation, obtained by 
using the rectangular integration technique, is the common 
positional form of the ideal PID algorithm which can be 


expressed as: 


where K, = Ke(T/Tr); ko = Kelip/T) 
The incremental form of the ideal PID algorithm results by 
taking the difference between the positional form of the 
@egorrunm at ‘the kthsamplle jnstant-(cf Equation? S622) and 
ene» positional form “for@*the’=(k-1)th sample instant. The 


incremental form of the ideal PID is then expressed as: 


u(k) - u(k-1) = (Ko + Ky + Kp)e(k) + (-Ke -2Ky)e(k-1) + 
Kye(k-2) (363) 


In the z-domain, the incremental form of the ideal PID algo- 





rithm is 
U(z) ane lan zc heta zs 
engel and’ lie oypelimay (Ox6e4) 
E(z) heer rani 


where ay = Kc + Ky + Kp 
ema e cat toh) 


a2 = Kp 
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4. Distillation Column Models 


4.1 Introduction 

Both linear and noniinear models have been developed to 
described the behaviour of the pilot scale distillation 
column used for the experimental phase of this work. Linear 
models have been used primarily in conjunction with the 
design of digital control algorithms while nonlinear models 
have been used to simulate the behaviour of the actual dis- 
tillation column for evaluating the performance of different 
icc alienahgornuthms. Wnepthemnextecsect ion, Section 4.2, the 
development of the nonlinear model of the binary distil- 
lation column used in this work is presented. In Section 
a3 Tse eGGet i Ong sot Bilsinvean | model sty forencontrolyad csit udsi esi gas 
discussed. Evaluation of distillation column models and the 
discussion of Simulation results are presented in 


Sectionse4.4,and 4.5,respectively. 


4.2 Nonlinear Model 
Thed \dynamicamodel,,;ofssthe hbinary distillation column 
used in this work was based on both the material and energy 
balance equations. Program development was directed to 
1) improving agreement between the simulated res- 
ponse and the experimental results; 
2) achieving flexibility of adapting different 


control schemes; and 
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3) improving efficiency in computational accuracy 
and speed. 

The original version of the nonlinear distillation column 
model used in this study was developed by Simonsmeier (50), 
engmesubsequent|y modified by Bijec (4). “In this work, 
modifications to improve the computational accuracy and 
Speed were introduced by both the author and his colleague 
Vee tr eusOnin #1978.4 Further jmodifiicattons tio, enhance,.the 
program's flexibility for studying the control behaviour of 
different control schemes and improve the agreement with 


experimental results were completed in 1980. 


4.2.1 Model Development 
The nonlinear binary distillation column model is based 
on the assumptions of 
1) perfect mixing in each stage; 
2) mass transfer and heat transfer equilibrium 
achieved instantaneously; 
3) negligible vapour mass holdup; and 
4) constant liquid mass holdup. 
From a macroscopic point of view, any physical system can be 
described by material and energy balance equations. This 
means that the basic balance equation 
Input + Generation = Output + Accumulation 42000 
must be satisfied at all times. Based on the aforementioned 
assumptions, the following sets of general equations 


describe the behaviour of any binary distillation column 
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Total Mass Balance 


Ltd Seem ne ntles eon coon 0 (4.2.2) 


Component Balance (on the more volatile component) 


dx, 


Hye 





een el <nt lee Cn ne enone nly nei) = onkn + 


he (4283) 


dt 


Energy Balance (on the more volatile component) 





dh 
ae — Senter — ena nn t Yn-18n-1 > Sn n * 
Fohe - Qin + Qin (4.2.4) 

where n = stage number 

L = liquid mass flow rate 

V = vapour mass flow rate 

S = side stream mass flow rate 

F = feed mass flow rate 

f = associated with feed 


WT = liquid mass holdup 

x = liquid composition 

y = vapour composition 
= liquid enthalpy 


H = vapour enthalpy 


Q, = heat loss 

qr = heat input 
A tray described by these equations iS shown in 
Figure 4.271% [he following functional™ relations are salso 
employed to simplify the problem formulation of the binary 


distillation column model 


Hy, = H(Yn) (42Ses) 
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Figure 4.2.1 ‘Schematic Diagram of an Equilibrium Stage 
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N= (Xp) (ames 6.) 
Hye Heer) (4.22. 17) 
Yn = Yher + En(Yn -¥Yn-1) (4.2.8) 
where y* = equilibrium vapour composition 

E = Murphree vapour efficiency 


For the pilot scale distillation column in the Department of 
Chemical Engineering at the University of Alberta, there is 
no side stream drawoff and the only feed to the column is at 
mieei uth stage (fourth, tray), so it tollows from Equations 
wae 2 tO 4e4-4 that the equations for tne nonlinear binary 
distillation column are 


AG tnemrebovier 








est = G (4.2.9) 
note: L,; is the bottom product flow. 

dx, 

dh, 


where Qe is the reboiler heat duty 
At the feed stage 
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QL5 


Sh 
C 2) 
4g = 
3 1S! 
f 

rr" 


2 i 
5 
- rt aba | 
' 
¢ 
~ 
4 
- a 
’ ee. ! 


' a ‘ 
frie 
tu ty 


a | 
- 


t*ow641dD megrts2 ebiz.an 


} “DOT : 


='S a 


feoe toltq ond tame 


resnten? [fsoltmen® 
_ 


iocss oft JA] ae 


(ax)¥ = x 

fu = 9% 
= “y s1edw 
= @7 









ooese Asti? ana 
,@ oF $. 5.8 
sesetinsent 


exb 
s TW 
= ae { 
au 
2 
: ™s 
th . 
ot 2d 979hW 
St a" 


oest sit \sA - 


+ ad - > 


43 


At the condenser 














Vg = Lig - Re = 0 C4821 5) 
note: Lig is the top product flow. 
xs 4 
WT 19 aa = Vo¥9 - Li10x10 - Rex1o0 (Ati 70)) 
dhig 
WTig Ta Ck un DL Oet Re tore UL 10, ~)Yeoo] ing 
(Ae eg) 
moc stages n= 2°3,4,56,/.8.9 
Lit] - Ln - Vn + Vn-1 = 0 (4.2.18) 
dx 
WT — 
ares SeLnriXnciese cnkne = \nynet neiyn-1 9 (452.19) 
dh, 
WT a 
n At = Ln+1Nn4]1 = Lyhn = VnHn te Vn-14n-1 (4 2820 } 


Note: for n = 9, hp,, the reflux enthalpy rather than 


hjo, the liquid enthalpy in the condenser is 


used. 

There are total of 30 equations available to describe 
the distillation columns dynamics, however, there are only 
29 basic unknowns, namely 10 liquid compositions, 10 liquid 
flow rates, and 9 vapour flow rates. Since the reflux flow 
rate is the manipulated variable for top composition control 
and reflux enthalpy is maintained slightly lower than the 
liquid enthalpy of the 9th stage to prevent flashing when 
the reflux flow is return to the column at the 9th stage, 
Equation 4.2.17 does not need to be utilized directly in the 


development of the nonlinear distillation column model. 
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4.2.2 Problem Formulation 

There are two approaches, namely the substitution 
method and the matrix inversion method, to solve the afore- 
mentioned system of equations. Previous work by Simonsmeier 
(50) has demonstrated that the matrix inversion method is 
more efficient than the substitution method, therefore, only 
the matrix inversion method will be discussed in this sub- 
Section. 

In the matrix inversion method, the system of equations 


for the distillation column is formulated into the following 


matrix expressions. 


ep 
uke (4.2.21) 
Le] 


x = BX (4. Omi. 2) 
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imtemdectatteds structure of Equations 4.2.21) and 4). 2 322 
is shown in Tables 4.2.1 and 4.2.2 respectively. From Equa- 


mona. 2.21, 1G TOlLOWS that: 


L “afr 
= hepyeiy [a ae (4.2.23) 
fg ee 


In Equation 4.2.22, x can be approximated by the backward 


i berence method. 


ee a (4.2.24) 


where k designates the kth integration step 
However, the right hand side of the above equation is a bet- 


ter approximation of Xavg Chat ex eeoU lee RUS, 


$h 


_ worgetusrot me fdor4d $.S.8 
16 svent 


} *- 5 ft > oi : soos sw 7 
norguss P]edue a i\ J /i 35 ,2anoaer ig 5 v 


- ve -_ esem off baer bossem 
cgrots. sat 6! 2 ih ‘wee alae r 


— 


rsrometvom!(2 yd . rmorveat 2rof pote to mez2¢2 nesnet Iast 
~ ; Ff ; -~tznomeb 2st (Ge) 
sorevears , bon nersusy Feu id non TsfotTTs sym! 

La : : ; Haenstag Ag'eéye xr1768 saz 


norjos¢ 


‘ 
a bs 
hag ‘ yu Da! f i aIaet ant 


Soe $[.5.# < d ‘+ nwode 2t 


~ 


thi rs ‘r .I9.8.32 were 


= 


; 
<-— = y 
Boeeenee Bhd yo veremixorgqs SO NMe9 A 55 . not tebp3 onl 


= 


«bon 29m s2natotttb 


nef 


t 





Tables ae 2.1 


SCUUGCURG OTrmrqudatiOne4 .2 1) 


ay -1 
=i 1 -1l 
=e al We eal 
el T=] 
=) }} Isat 
aj at 
stati 1 -1l 
“fh 4] ] 
=1oeet 
-1 
a) by -1 
a2ub2 C21 
23 b3 C3 =i 
oa ba Rae: 
#5 bs Cham =i 
&—6 bE ce =k 
a b7 Cae ete:l 
4g bg cg 
#9 bg 
where a; = -h, / H; i = 
By a tiie iy le 
Crikey steer ee eat bea 
dh 
ope A + Cine Op) mate 








(=) {eh Key =) 


‘~ 
a 2 


+9 S°1NT30'. z 


™ — « 
» TGissae) 





—— + 





—— tr nae 





Tabdes 4s. 2292 


SUGUCCUT CMO Tm LOUd COM eo sci. 2 





oh cpesed 
bo a9 Co 
0 Boag 
0 0 ba 
0 0 0 
=—o-——9* Co 
0 0 0 
0 0 0 
0 0 0 
he 0 0 


=) 
iT} 





a5 


8 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 Q 0 0 
ee wl 0 0 
Wied fe 0k! 0 


b6 


Doyo Os Gs 05, LU 


i 


0 be ag Cg 
OMM COMED G nao 


LVy (1 - Kn) 2 Ln+] os Vn-1J 








us 
0 m1 Te 
oe 
0 X3 
0 X4 
0 X65 
0 1X6 
0 X7 
0 XQ 
Cg X9 
410 a 


_ [Vs (ey - Ks) - nas - Va + F (R - a 


WT5 


1 
WT Wire ene 





1 





Ln+] 
WT 


Yn / Xn 


LG 


nS? 2E3y AEs Orta sts. cU 


n 


V2, SFA SSr, OL ESe9 


46 


os 





: ' > 
5- Se" 8 re 
. 1 40 940 Pouti2 
S.S.h a suoa to avurs 


$79 


to 


™ 


ex f ' n q 0 f) f} S 2 7 &] 
|. #4 | , 
f) : 7 } , rf cf F 
_ ies } 4 ] J u * U Vv - S & ¢* 
| - . : . 
“ 0 G 6 0 6) 0 63 re ee a) | 
pe | 6 0 Oo so pe ope Q 6 












3 Pew) 4+ gv = ad - fee OD at) 
i | ’ 2TW 


[ 
ae 8,V,5.¢.4 £.5 a a) fant yen’ = 
a mig . 
i > - ad >. 
i 











47 


; avg (x 


" 
S 
o1 
2) 
ca 
Ww 
|>< 
CR Oo 
+ 
w 
oe 
I 
x< 


If the integration interval is several orders of magnitude 
smaller than the mean time constant of the process, the fol- 
lowing expression will be valid. 

BR—~ BK-1 (40226) 


Therefore, 


= 0.5(BkK-1xk + gk-1xk-1) Gamer 7)) 


eee CIPS HPS CEBKAN)- Min +siongatBhel yk wh (4.2.28) 


4.2.3 Numerical Solution 





There are numerous computer library subroutines for 
solving system of equations using direct matrix inversion, 
however, these subroutines are intended for general matric- 
es. For matrices with a large number of zero elements, such 
as the tridiagonal matrices, direct matrix inversion is not 
necessary. For inverting tridiagonal matrices, the Thomas 
algorithm (24) is the most effective and efficient approach 
and-LittHs chosen. for solving thesiiiquidsjcompositions@inethe 
model. For matrices which have a large number of zero ele- 
ments but do not have the structure of a tridiagonal matrix, 
an iterative approach could be more efficient than direct 
inversion method if the matrices are stable and good initial 


guesses are available. The iterative method is particularly 
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Suitable for distillation modelling, because the distilla- 
EAOneCOLUMNR ALS baeS WOW DHOCES Suastinusasthe,nnformationy, from the 
previous integration cycle will be a good initial guess for 
the present integration cycle, if the integration interval 
is small enough. Among iterative methods, the Gauss-Seidel 
method is very effective so it is chosen for solving for the 
liquid and vapour flow rates in the model. 

USing athena Thomas sal,gorithm gto sisodves) for) these liquid 
compositions and the Gauss-Seidel method to solve for the 
liquid and vapour flow rates, the computational efficiency 
has been improved drastically while the accuracy of the 
simulated results remains the same or even better. For the 
solution procedure employed by Bilec (4), the amount of CPU 
time required by the Amdahl 470/V7 computer for a typical 
Open loop simulation (total of 360 inteyration steps) was 
more than 9 seconds while the new solution scheme used in 
This work requires yless than.3 seconds; of.CPU,time, for the 
Same simulated time period. The FORTRAN program listing of 
the new model and the associated subroutines are available 


from the DACS Centre. 


4.3 Linear Model 

Nonlinear models are very useful for simulation studies 
of physical systems, however, to use nonlinear models, such 
as the nonlinear binary distiliation column model, in con- 
trol algorithm design is inappropriate simply because mathe- 


matical manipulation is almost impossible. Linear models or 
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transfer function models, on the other hand, are favoured by 
control engineers for their mathematical simplicity as well 
as for the ease of determination of model parameters. 

The most commonly used transfer function models in 
classical control system design are in the Laplace domain, 
because they represent continuous systems. Most industrial 
Processes can be represented by first order plus time delay 
models or second order plus time delay models. Higher order 
models, are seldom used because the parameters in these 
movers) are dttticult toltidentify bye simple eprocedures. In 
this work, only first order plus time delay models are con- 
Sidered because they simplify the development of digital 
control algorithms. 


—_—— 


4.4 Evaluation of Distillation Column Models 

The use of nonlinear model simulation is to provide 
control engineers with a more realistic evaluation of the 
performance of the control algorithm prior to the actual 
implementation, therefore the simulated open loop responses 
Should match closely with the .actual experimental data. 
Three sets of open loop tests were conducted on the pilot 
scale distillation column for verifying the nonlinear model. 
The tests were performed for +25% step disturbances in feed 
flow rate, +10% step disturbances in reflux flow and +#8% 
Step disturbances in steam flow rate. In order to match the 


Simulated open loop responses with the experimental results, 


the following procedure was followed: 
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Obtain all variables of the nonlinear binary 
distillation column model at the initial steady 
State acondneians:: 

a. Feed flow rate, reflux flow rate, and steam 
flow rate. 

b. Feed composition, top product composition 
and bottom product composition. 

c. The column temperature profile. 

d. The liquid composition profile. 

é€. the *liquid mass holdup profile. 

f. The temperature, liquid composition~= and 
liquid mass holdup profiles are not critical 
to the simulation if only the terminal com- 
position responses are of particular inter- 
est. However, the availability of actual 
data for these profiles will ease the veri- 
ficabicne pyocedure. Pinpotnis, project, data 
from Bilec (4) are used as the initial esti- 
mate of these profiles. 

Estimate’ the »initial eheat’ loss ) profiles p0ne 
method °of ‘estimating tnisSproftidteryisetoetdtstri= 
bute the total heat loss (obtained from the 
experimental heat loss) to each stage according 
to the difference between the stage temperature 
and the ambient temperature. 


Estimate the initial tray efficiency profile. 
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Enter the above data into the appropriate loca- 
tdonsanéthe datiat Tile. 

Simulate the steady state operation of the dis- 
tillation column using the above data. 

If the simulated steady state conditions, i.e. 
the top and bottom product compositions (since 
the terminal compositions are of primary inter- 
est), do not shown close agreement with the 
actual experimental results, the heat loss pro- 
file is adjusted and Step 5 is repeated. 

If the steady state conditions are in agree- 
ment, then the simulated conditions will be 
used as: ithemianitral @conditions foreitheropen 
loop simulation. 

Estimate the above profiles for the new steady 
Statviesycond i trons. 

Simulate the open loop response. 

If the simulated new steady state conditions do 
not show agreement with the experimental steady 
State conditions then adjust the heat loss pro- 
file for the predicted new steady state condi- 
tions*andsrepeat Step. 

If the transient response of the terminal 
compositions are not consistent with’ the 
experimental response then adjust both the 
initial and final tray efficiency profile and 


repeat Step 8. 
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12. The model parameters must continue to _ be 
adjusted so that the steady state predictions 
agree with the experimental data before match- 
ing the transient responses. 

Since bottom composition control study is the primary 
metcern Ot this project, the matching of top composition 
responses is not of concern. The simulated nonlinear model 
predicted bottom composition responses using the above pro- 
cedure and the experimental responses are shown in Figures 
4.4.1 to 4.4.6. The simulated responses shown reasonable 
agreement with the actual experimental data except that the 
nonlinear model predicts faster column dynamics than is 
Observed experimentally when the disturbance is either the 
Steam or feed flow rate. 

In order to perform simulation studies using the non- 
Pinear distillation column model, a set of transfer func- 
tions that can represent the simulated nonlinear’ model 
response are required for control algorithm design. Since 
the prime objective of this work is the study of the per- 
formance of various control algorithms for bottom composi- 
tion control for the column subjected to a feed flow rate 
disturbance, the transfer functions of interest are those 
relating bottom composition and steam flow rate and between 
bottom composition and feed flow rates. The simulated non- 
linear model responses were approximated by first order plus 
time delay transfer function models. Selection of the trans- 


fer function parameters was performed on a trial and error 
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Figure 4.4.1 Comparison of Experimental and Pre- 
dicted Nonlinear Model Responses of 
che=eBOtrom “COmpOoSTtion “TOT "a= 25% 
Step Decrease in Feed Flow Rate 
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Figure 4.4.2 Comparison of Experimental and Pre- 
‘dicted Nonlinear Model Responses of 
thes Bottom «eComposmition Fifor fate 25% 
Step Increase in Feed Flow Rate 
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basis with the aid of a small IBM 360 CSMP simulation pro- 
gram. The transfer function models are listed in Table 4.4.1 
and the comparison of the linear and nonlinear simulated 
bottom composition responses are shown in Figures 4.4.7 to 
4.4.10. These figures show that excellent agreement between 
the linear and nonlinear model simulated bottom composition 
responses was achieved and the linear models, SLMO1 to SLMO4 
can be assumed to have minimal model parameter error in 
representing the noniinear simulated bottom composition 
response. 

Similarly, the experimental open loop results for the 
+25% step disturbance in feed flow rate and +8% step distur- 
bance in steam flow rate were fit by first order plus time 
delay transfer function models for use with the control 
algorithms to establish the control law for use in the sub- 
sequent experimental evaluation of the performance of the 
different control algorithms. These transfer function models 
are listed in Table 4.4.2. Excellent agreement was also 
achieved between the experimental and simulated linear model 
responses as shown in Figures 4.4.11 to 4.4.14, hence models 
ELMO1l to ELMO4 are assumed to have minimal model parameter 
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SLMQ1 
SLMO2 
SLMO3 
SLMO4 


Table 4.4.1 


Transfer Functions for Simulated Open 
Loop Responses from Nonlinear Model 











Input X(s)/U(S) Ky (wt.%/g/s) Tp (s) 
-25% FE XB(s)/FE(s) 1.00 246 
+25% FE XB(s)/FE(s) 2.26 577 
-8% ST ee ene) bree 4 563 
+8% ST XBls)/ST(s 4829 633 


* - measurement time delay is not included 


Model 





ELMO1 


ELMO2 
ELMO3 


ELMO4 


Table’ 4.4.2 


Transfer Functions for Experimental 
Open Loop Responses 





Input Mis) LUGS) Ky (wt.%/9g/s) hy 13 
25% FE XBNOS)/ FES) Tea 0 415 
a rAoy apes ese tes 2.20 929 
=9% Slee XB(S)/SI(S US wed (a) POSS 
TO Sie AX US) /-S1 CS) -4.29 918 


* - measurement time delay is not included 
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Figure 4.4.8 
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Comparison of Simulated Nonlinear 
and Linear Model Predicted Responses 
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Figure 4.4.9 Comparison of Simulated Nonlinear 
and Linear Model Predicted Responses 
Oe LUCE bOLt Om .COmn0s. ion, for a 8% 
Step Decrease in Steam Flow Rate 
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Figure 4.4.11 Comparison of Experimental and Pre- 
dicted Linear Model Responses of 
Bottom Composition for a 25% Step 
Decrease in Feed Flow Rate 
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Figure 4.4.12 Comparison of Experimental and Pre- 
dicted Linear Model Responses of 
Bottom Composition for a 25% Step 
Increase in Feed Flow Rate 
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Figure 4.4.13 Comparison of Experimental and Pre- 
dicted Linear Model Responses of 
BOCCOMMmECOMDOSTIULLOM mtOb sd) Gp “Step 
Decrease in Steam Flow Rate 
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Figure 4.4.14 Comparison of Experimental and Pre- 
dicted Linear Model Responses. of 
Bottom Composition for a 8% Step 
Increase in Steam Flow Rate 





5- Digital Simulation Studies of Bottom Composition Control 


5.1 Introduction 

The objectives of studying bottom composition control 
using on the nonlinear distiliation column model are 

pros dain an insight Into tne control behaviour (of 
the ®distillation column under different control 
algorithms. 

2) To compare the performance of the different con- 
trol algorithms under both servo and regulatory 
control. 

Sym lO Sgaini lan cinsight#into ‘thew abidity vof ithe tdifi- 
ferent control algorithms to compensate for pro- 
cess and measurement time delays. 

4) To gain experience on the tuning techniques for 
tne wWatdenentvicontrol algorithms: 

The schematic diagram of the distillation column unit 
to be simulated is shown in Figure 6 2.1 (cf. Chapter 6) and 
umes description tof étherisystem isitgivenrineSection 642 lof 
Chapter *6." The digital’ simulation: of ‘the “control “system was 
carried out on the Amdhal 470/V7 computer using the non- 
linear distillation model described in Chapter 4. Bottom 
composition control behaviour was simulated for step changes 
Ofe 2252 cin feed flow rate and +1% changes in composition set 
point. The sampling time was 3 minutes consistent with the 
measurement delay of 3 minutes. Simulations were performed 


for a total of 165 minutes with the load disturbance or set 
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point change introduced 15 minutes from the start of the 
Simulation. These conditions were also used in the experi- 
mental evaluation, so the simulated results trom this sec- 
tion can be compared directly to the results from the exper- 


imental evaluation. 


5.2 Performance Criteria 








In most process control studies, error-integrals are 
used as the performance criteria for evaluating control 
algorithms. These measures are used because they character- 
ize the entire time response of the controlled variable when 
the system is subjected to any disturbance, unlike other 
Criteria, such as the percent overshoot; decay ratio; rise 
time, response time, etc , which are useful for evaluating 
the system performance for servo operation. The value of in- 
tegral of the absolute value of the error (IAE) was used as 
the criterion for evaluating the performance of the differ- 
ent control algorithms in this work. Although the perform- 
ance* of a’control algorithm'is judged ‘primarily by the IAE 
value calculated on the response of the controlled variable, 
other practical considerations such as the response of the 
manipulated variable, e.g. valve constraints, maximum rate 


of change, should not be ignored in actual implementation. 


9-3 Performance of the Algorithms for Servo Control 








Since the distillation column is a highly nonlinear 


process, especially the response of the bottom composition 
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to a change in feed flow rate, it is necessary to use diffe- 
Ventesets Of Comenol er cangtants! forucontrolling different 
gusturpances. Thetinitial controllerSsettings for the PI and 
PID algorithms for servo operation were calculated using the 
correlations suggested by Smith (51) based on either the 
BemO03 or SLMO4. model-s).(cof. Séctiion -4.4).-The -IAE?valués?ror 
the performance of these control alyorithms, after tuning to 
minimize the IAE value, are summarized in Tables 5.3.1 and 
5.3.2. The bottom composition control responses using these 
@eecmrcnms are snown in- Figures 5.3.1 to 5.3.4. 

For the deadbeat algorithm, the Dahlin algorithm and 
the Smith predictor, the SLMO3 and SLMO4 models were used in 
the design of these algorithms. The IAE values for the per- 
formance of these algorithms are summarized in Tables 5.3.1 
and 5.3.2. The control performance achieved using these dif- 
ferent algorithms is shown by the responses given in Figures 


Dea tO 5.3.40. 


9-4 Performance of the Algorithms for Regulatory Contro 





Bor the#Pl and P4D contol algoritims23the? initliat, sets 
oF seontroller-cConstantsosare; again, cabculatedS.using 7the 
correlations from Smith (51). The IAE values that resulted 
Using theiPlwand.PID algorithms, fom regulatory control are 
ha S tedé aim) Tabl espe5. 4. 1° and “S¢4. 20 ewni ie the - responses, of 
bottom composition using these two control algorithms are 
presented in Figures 5.4.1 to 5.4.4. For the deadbeat, Dah- 


Tin and Smith predictor algorithms, the initial controller 
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Tabress:o7 | 


Summary of Simulated Results for 
Controlzor bottom Composition for 
a@ rie Step Change in Set Point 








Algorithm ne KI KD Kp Tp Td IAE Figure 
PI -0.160 -0.039 — — — — 966 5.3.1 
PID -0.198 -0.045 -0.192 — — — 758 5.3 3 
Deadbeat — — — OC OS Ol mee) cae Oe 
Dahlin A = 14s — —_ Seo A DO Om O? o).2e ad 37, 
Smith Pred. -0.678 -0.210 — Biyhyaih Rise) aks “7h Gn sia!) 
Ke = (g/s)/wt.% Kp = wt.%/(g/s) 
irr (g/s)/(wt.%-s) Tp = Ss 
Ky = (g/s)/(wt.%/s) Tq = 8 
IAE = wt.%-s 
Table ws. 3.2 
Summary of Simulated Results for 
Control of Bottom Composition for 
a -1% Step Change in Set Point 

Algorithm KG KI KD Kp Tp Td AE Figure 
Pl -0.160 -0.039 — — — — 949 5.3.2 
PID -0.240 -0.053 -0.197 — — — 740 5.3.4 
Deadbeat _ — — =4°29 6339 137 (67. 2aeon S60 
Dahlin » = 20s — — =455 2921633010, / Od me OO 
Smith Pred. -0.948 -0.263 — aA C9 053 5 Oy Osan re 
Ko = (g/s)/wt.% ay wt.%/(g/s) 
KT @= (9/5) /(wt.%-s) Tp =s 
Ried s)/ (WE e as) Td = 8 
IAE = wt.%-s 
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Figure 53. PR - Algorithm, Control) of Bottom 
Composition for a +1% Step Change in 
Set Point (Run S-PIO1; he = <=02160; 
Ky = -0.039) 


“wo 






= : ee —— 
. 
- vw Ar eLe., 
4 ol . 4 ~ S 
a - e r 
+ — 
, i€ a4 » 4 
OIA NT ODIO DIODE aly 


{> 


Oro RSD 





14.50 15.50 


ole G/S9 


13-50 


£2260 





HUE Eee a) 10) - 


8.00 


(WIZ) 


“H'00 


0.00 40.00 - 60.00 , 
TSE E 3S?) x10 


Fugure 322 "Pe Algorithm  €ontral O ia OL On 
Composition for a -1% Step Change in 
Seteeporit Ss RUNT o=Pl0z: Ke = -0.160, 


Kr = -0.039) 








67 


-00 20.00 40.00 60.00 80.00 100.00 


80.00 100.00 





Gomeneraneceie xe Do seSIOIOU UD 
‘, * : 


15.50 


(G/3) 


T 


UD 


(WT) 


14.50 


-50 


™m 
4 


oO 

a 

ra) 

b.00 20.00 40.00 60.00 , 80-00 100.00 

HIME =ESEG) “xeho 

oe 

=} 

ioe) 

= 

eS 
© 
© 

i) 

CO 

+ 


Co | 
= , 
“).00 20.00 40.00 «60.00 80.00 100.00 


PUNE eSEGs x 107 


Piqure 5.3.3 PIO Algorithm .Control of Bottom 
Composition for a +1% Step Change in 
Set Point (Run S-PIDO1; Ke = -9 198, 


Keer. 045, Kp 4 +02192) 








68 





(G75) 


oT 


14.50 15.50 


13.50 


69 





i 

w . 

b.00 0.00 40.00 60.00 , 80.00 100.00 
THEME testCy) *10 

F | 

oS: 

wo 









0.00 40.00 60.00 , 80.00 100.00 
FRE (68SEC) x10 


Pui gtres 54344 np Pil Die Al gonmithmbepControl; of » Bottom 
Composition for a -1% Step Change in 
Set Point (Run S-PID02; Ke = 0.240, 


K; = -0.053, Kp = -0.197) 


-00 


» CF wud ae eer oo) 
ae “> > . . 
7’ > 7 
- : 
7 Ve ' 4 
-_ 


od caee as ~ 





ST (G/5) | 
15.50 


12-50 


14.50 


13.50 


70 





-00 20.00 40.00 60.00 80.00 100.00 


TIMES SBC)) ql 0) 


Qo DPODOOGODHOOOHO OOO DODO DOOD OOOO OHO DIO 





-00 


F tquce5 2325 





40.00 ~60.00 80.00 100.00 


0.00 s 
Tal Mee DeSitatee) x10 


Deadbeat Algorithm Control of Bottom 
COMPOS LG uOtel Ored tle sep Change. in 
Sete Paint. (Run S=DBOL; Kp ae 5a te 
Ty = 0563205 Tq = .187.0) 








ue me 
ho 4 ‘ 
- 0d. Ef 


| 
(4 P| 
\ 





yr 
t 


ang 


idl 


15.90 








O 

uD 

N 

™.00 20.00 40.00 60.00 , 80.00 
ecomomy) x10 


100.00 





© 
S 
amate: 20 Abe0er 4 60.00 80.00 100.00 


-O0 ; 
cl ViESastee Ese 8 x10 


Figure 5.3.6 Deadbeat Algorithm Control of Bottom 
Composition for -1% Step Change in 
Sec GhOEnt ((CRUN ws-Ub02) GKR = 4°27 OF 


= P 
Ty = 633.0, Ty = 187.0) 


y {2 





02.61 


os 








0) 


ODOR 


Lab 8 43 


CRC LONOEX 


¥ Ty 


15.50 


7 a) 
14.50 


13.50 


om 
8.00 ol2-50 


6.00 


(WT) 


° 
mo 
<2 


2-00 





-00 20.00 








Beye 


-00 20.00 


TIME 


Pigure: 5.3. 7 


Dahlin 





40.00 60.00 


(SEC) *10° 





40.00 - 60.00 


(SEC) xo) 


Algorithm Control 


12 


80.00 100.00 





80-00 100.00 


Of “Botton 


Composition for +1% Step Change in 
Set Point (Run S-DAHO1. 


Tet. 
uk 


Ke =D. 2e, 


DomUs dq = leg.0s a= 14.09 








78 


15-50 


14.50 


351. GAo4 


13.50 






Cy 
Up) 
w 
™b.00 20.00 40.00 60.00 , 80-00 100.00 
PAE NERC re 
5 
== 
00 


© 


OD NOOO OT ID IIS ISIS IIIS III ISIS IIIS ISIS IIIS 


2) 

‘2 

“b.00 20.00 49200R) 60400 , 80.00 100.00 
rep) Comat e circa: x10 


Figure 5.3.8: Danian Al\gorithm * Control of awbottom 
Composition for -1% Step Change’ in 


Set Point (Run S-DBO2- Ky = -4.29, 
Tp eouen0n) 1G = 187.0, x = 20.0) 





15.50 


14.50 


pier 
-50 


gy™ 


12-50 


=) 
© 
ee) 





00 20.00 
T 


-00 20.00 
ir 


Fitguire 57.3.9 


itv 


IME 


Smith 
Compo 


Set Pome 


Te. = 
Pp 
KT = 





PSEC 


40.00 . 60 
CEG.) 


Predictor 


74 


40.00 60.00 80.00 100.00 


wis 





-O0 5 80.00 100.00 
*10 
Gomluro.- “OT "Botton 


STCTOn Om ie stotep, change. im 


563% 0% Tq = 
-0.210) 


(Run stckr Olea ke © =c1oRy 


187.0, ke = -0.678, 





(WTZ) 


15.50 





-O0 20.00 





SIME 


40.00 
UOBE) 





Composition for 


SeiteeP.0 1 it 
LP i Css. 0; Tq = 


I 


-0.263) 





75 





60.00 


“See 


pL ME wip ese *10 
Figure’ 5.3.10 Smith Predictor ‘Control 


80.00 100.00 





80.00 100.00 


or. Botton 


Sie ocep eChange cin 


(Run S-SLPO2,; 


way A'F29% 


ioe 0 Ke 07 945) 


ra ee eae 
* al Py 


— . 





nATOOCKy 


oe 
— 





eres 





aby e. 5. 4:4) 


Summary of Simulated Results for 
COnL GO or Scottom Composition for 
a -25% Step Change in Feed Flow Rate 


Algorithm x KI Kp Kp Tp Tq 
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Summary of Simulated Results for 
Control of Bottom Composition for 
a +25% Step Change in Feed Flow Rate 
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parameters used for the servo control tests were those used 
fomgecne regulatory control tests. Since unsatisfactory per- 
formance resulted using these parameters, an improved tuning 
procedure especially for these algorithms was developed. The 
performance with each of these algorithms as indicated by 
the lower IAE values, before and after the new tuning proce- 
dure was employed, are also given in Tables 5.4.1 and 5.4.2. 
The controlled composition responses using these algorithms, 
with and without tne improved parameters, are shown in Fig- 


igecmeae4. 5) tO 5.4.16. 


9-5 Tuning Techniques 


Tuning techniques tor conventional PI and PID control 
algorithms have been reported by numerous workers (cf. Chap- 
Gecere). in this work, the correlations reported by Smith 
(51) were used to provided the initial sets of controller 
GComstants for both the PI and PID control algorithms. The 
procedure for adjusting the values of the constants to allow 
for the effect of sampling, as described by Smith (51), was 
followed. The final tuning constants were obtained by trial 
and error using the IAE value as the guide. 

In the case of the Smith predictor, the initial values 
of Ke and Ky were calculated using the same models as used 
forpthe controtji.enm constantaicalcubatwoan Tor SthesiPl. control- 
ler except that the system time delay term was deleted. 


Again, the correlations from Smith (51) were used and the 
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tuning procedure followed was that of minimizing the IAE 
value. 

POGmcnemsDanhin= algoripthm, the initval tuning involved 
Setting A, the time constant of the desired closed loop ser- 
vo controlled response, to be the same value as the value of 
Th, the time constant of the process. In the subsequent tun- 
ing, the desired performance is obtained by adyusting \ with 
the control performance improving as the value of A is 
decreased. The eftect of varying the value of A on the 
response of the controlled and manipulated variable is 
MeMlce rated in Figure 5.5.1. 

From the digital simulation studies, it had been found 
that conventional tuning techniques as described above can- 
not be used to further tune the parameters of the Dahlin 
@ucgrmunm and the Smith predictor’ to yield satisfactory per- 
Mocmance tor regulatory control. This is shown by the IAE 
values in Tables 5.4.1 and 5.4.2. Since the deadbeat algo- 
rithm has the same general structure as the other two algo- 
rithms except that no provision exists for convenient adjus- 
ting parameter values, the performance of the deadbeat algo- 
rithm is governed strictly by the process model parameter 
values used in the algorithm. A sensitivity analysis on the 
influence of model parameter values on the performance of 
the deadbeat algorithm, can provide guidance as to the type 
of tuning procedure that might be appropriate for improving 
the performance of all three algorithms. Since the time de- 


lay of the process can be estimated quite accurately, only 
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the influence of the values of the process gain and time 
constant on the performance of the deadbeat algorithm were 
Studied. A summary of the results from the sensitivity anal- 
Meseson the deadbeat algorithm is given in Table 5.5.1. As 
Shown by these results, the model parameters used in the 
deadbeat algorithm can in fact be treated as tuning parame- 
ters. The performance of the deadbeat algorithm improves if: 

Hyethe process gain; iS gunderestimated, upg to’ 15% 
lower than the true value; and/or 

2) the process time constant is overestimated, up to 
15% higher than the true value. 

On the basis of these observations, the following tun- 
ing procedure tor regulatory control is recommended: 

1) adjust the tuning parameters in the algorithm, 
esd. ere in whelwWah lpn raircgorithm Sruntiiyl inoofurther 
improvement (for the deadbeat algorithm, this 
step will be omitted); 

2) adjust the process gain constant used in the 
model to improve the response of the system until 
no further improvement is possible; 

S)sadjust, the process time: constant ‘used “ing =the 
model in the same fashion described in Step (2); 
and 

4) repeat steps (1), (2) and (3) if necessary. 

By adopting this tuning proceuure, it was found that by 


adjusting Ky and Tp of the process model used in the dead- 


beat, Dahlin and Smith predictor algorithms, a significant 
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Table 5.521 
Sensitivity Analysis of the Influence 


of Model Parameters on the Performance 
of the Deadbeat Control Algorithm (Ty = 187s) 


% Deviation From 
Actual Value =30%e 43-15% 40% +15% +30% 
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improvement in their control performance for regulatory con- 
trol resulted (cf. Figure 5.5.2). A summary showing the pro- 
gression of tuning the model parameters in the deadbeat 
muce@nemneusing this technique 1s “Given Tn Tabele’™. 522. A 
comparison of the transient responses of the composition 
control achieved using ditterent model parameters 1S shown 
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oo Discussion of Simulated Control@Ppertormance 

For servo control operation, using the IAE value as the 
performance criterion, the deadbeat algorithm gave the best 
control performance and the Dahlin algorithm provided almost 


the equivalent performance since it was tuned in a manner 
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Summary of Simulated Results for Deadbeat 
Rugonvent=contro lof bottom Gomposition for 
a -25% Step Change in Feed Flow Rate (Ty = 187s) 


RUN Ky (wt.%/g/s) Tp (s) IAE (wt%-s) 
S-DBO4 -5.24 sean. 5382 
S-DB30 -4.54 563 4594 
S-DB31 -3.78 563 3913 
S-DB32 aSe2u 563 3328 
S-DB33 =e 563 2832 
S-DB34 ONS? 563 2426 
S-DB35 =heoy 563 3869 
S-DB36 -2.44 563 2532 
S-DB37 32320 563 2480 
S-DB38 ene 591 2516 
S-DB39 22632 535 2408 
S-DB40 = 263 2 508 2408 
S-DB41 E2632 483 2408 
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Similar to that used for the deadbeat alyorithm, i.e. a very 
Small value of . The IAE value for Smith predictor control 
is very close to those of the deabeat and Dahlin algorithms, 
however, the response of the manipulated variable under 
Smith predictor control is more oscillatory than was the 
case for the other two algorithms. The deterioration in the 
performance of both the PI and PID algorithms is mainly due 
to the significant amount of time delay in the system (pro- 
cess and sample time of the gas chromatograph). However, it 
must tbe ‘noted, that derivatavesractiom provides «significant 
compensation for the time delay as shown by the much lower 
IAE value using the PID alyorithm compared to the value that 
resulted using the PI algorithm. 

In the case of regulatory control, using a conventional 
tuning technique, the performance of both the deadbeat and 
Dahlin alyorithm was interior to that of a tuned PID algori- 
thm when the system was subjected to a +25% step disturbance 
in feed flow rate. Furthermore, their performance was even 
worse than that obtained using the PI algorithm when the 
disturbance was a 25% decrease in feed flow rate. The per- 
formance of the Smith predictor for regulatory control oper- 
ation also deteriorated in a similar manner. Meyer (37) ob- 
served that the performance of the Smith predictor for regu- 
latory control is dependent on the ratio of the time con- 
stants of the process and load disturbance transfer func- 
tions even though the algorithm is designed to cancel the 


time delay term in the closed loop system characteristic 
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SGila vn sen UN SaeWOCKRIMIt ‘waseefoundmathat. forguregu latory 
control, the performance of the deadbeat, Dahlin and Smith 
predictor algorithms using conventional tuning techniques is 
unsatisfactory when the time constant of the process is much 
Kaggeristhanathat)of theyload transfers function:*Inetherother 
wonasr, upachecoked limofathesefiects ofs there odds ditstunbance fon 
the process controlied variables results beiore the control 
action takes place. Using the improved tuning technique for 
the deadbeat, Dahlin and the Smith predictor algorithms, 
this difficulty has been overcome since the improved tuning 
procedure is essentially to force the controller to calcu- 
Paces anStGOnGer. andmore rapid. control action by using a 
transfer function model, in the design procedure, which is 
less sensitive and has slower dynamics than the actual pro- 
cess. 

The simulation results showed that the response of the 
controlled variable (bottom composition) in the tests of the 
deadbeat, Dahlin and Smith predictor algorithms for both 
regulatory (after the use of the improved tuning procedure) 
and servo control were quite smooth. However, the response 
of the manipulated variable (steam flow) using the deadbeat 
algorithm might not be acceptabie tor industrial application 
(cf. Figure 5.3.6, where the steam flow exhibits more than 
100% overshoot). For the Dahlin algorithm and the Smith pre- 
dictor, this problem can be avoid by adjusting the tuning 
constants in these algorithms to obtain an acceptable re- 


sponse of the manipulated variable without causing a serious 
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deterionati onipun’tthettcontroll epenfionmanc eaatc fv ePigures: 5.4.13 
and 5.4. 4)" 

Although the response of the controlled and manipulated 
variables are quite sluggish and oscillatory when composi- 
Grom tsMcontrollede with ‘the? Pilo@andePplDealgorithmss for time 
delay processes, these algorithms do have the advantage that 
they are model independent and so are easier to implement 


for any control system. 












ts un 
a 
} @ 
—-<¢ 
! ‘e fortneo sity mi ae rise tysteb 
tc PCy . ane | } ° J . “ 4 is ee , ae 
: | (hI 4.2 Das i 
lovtng> af? Yor sena@es" Bae dpvonys tA Asa 
34 iil | | | . 
opooufe asétup sve geldstvev oF 
j ang nat bulflortnegs’ ef 
‘ pls Seensy ,~e¢3eeo9 " Y 
¢ rubtaqabatr (epom ots  ¢ 
; 


ie 
wageve Torres “ns Be 
; Re 


" 





6. Experimental Implementation and Evaluation 


Ot thiunt roduction 

The ob,ectives of the experimental work of this project 
were somewhat different from those of the digital simulation 
study outlined in Chapter 5. In the experimental work, 
efforts were directed to 

1) implementation of the different control algo- 
rithms using the process control software system 
Dal S01 (OR 6 263564165); 

2) evaluation of the control performance of the dif- 
ferent algorithms for a -25% step change in feed 
flow rate; and 

3) on-line tuning techniques for the different con- 
trol algorithms, particularly the tuning tech- 
nique tor the Dahlin algorithm. 

In the evaluation of the performance of the different 
control algorithms, the deadbeat algorithm was not consider- 
ed, primarily because the algorithm did not have any provi- 
Sion for tuning parameters to adjust different operating 
conditions iGcfrersection 394) * 

In théoinexta section, eS ecttiontes2A thes description eof 
the’ pilot scalettdistiliation ‘column *untt Pised® in’ this*pro= 
ject is presented. An outline of the basic operation of the 
DISCO system is given in Section 6.3. The control perfor- 
mance results for the different algorithms under regulatory 


and servo control operation are presented in Sections 6.4 
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and 6.5 respectively. The on-line tuning techniques used in 
the experimental evaluation of the different control algo- 
rithms are given in Section 6.6 and the discussion of expe- 


rimental results is presented in Section 6.7. 


Gees Description of the Distillatgvon=eolumn Unit 





The experimental evaluation of the different control 
algorithms was conducted on the pilot scale binary distil- 
lation column located in the Department of Chemical Engi- 
neering. The column originally designed and fabricated by 
Syreek (55) has been uSed/ for .dvanious «types Sof_/control 
Stuqies during theslastmfew years. The column jis) 22.9); cm in 
diameter and has eight bubble cap trays on a 30.5 cm spac- 
ing. Each tray contains four bubble caps. The column is also 
equipped with a total condenser and a thermosyphon reboiler. 
The liquid levels in the condenser and the reboiler are con- 
trolled by manipulating the top and bottom product flows 
respectively using local analog PI controllers. The column 
is maintained at approximately atmospheric’ pressure. by 
adjusting the cooling water flow through the condenser. The 
feed flow is pumped to the column at the fourth tray after 
passing through a preheater for maintaining constant feed 
enthalpy. Similarly, the-reflux is also preheated to main- 
tain constant reflux enthalpy before it is returned to the 
column at the eighth tray. A schematic diagram of the dis- 
Plstd ce PODMBCOLUI Udita Sao NOWN New hl GUGeO..2.1 and) the 


Lypicalm operating conditions used) during the experimental 
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muUNnSse anetiastedurin! Labihew6.2aliy Specificaedetaidsymnegarding 
the design of the binary distillation column unit are given 


by Sivrecek? (55)% 


aD em6 recan. 


Steady State Operating Conditions of 
Chee bia RyesD sic lel ata onerco umn aU nat 





Feed Flow Rate io..0uG/S 
Bottom Product Flow Rate Oo Ug ss 
Top Product Flow Rate G0eG /aS 
Reflux Flow Rate OR 0erG 7S 
Steam Flow Rate han Omdy/ S 
Feed Composition 50.4 wt.% MeOH 
Top Composition 95.3 wt.% MeOH 
Bottom Composition 5.0 wt.% MeOH 


In order to provide data acquisition and implementation 
of the different control schemes, the distillation column is 
interfaced with the distributed HP 1000 computer system. A 
HP 5720 gas chromatograph which is also interfaced with the 
HP 1000 computer system, is used for direct measurement of 
the bottom product composition. During each experimental 
run, eight key variables, namely feed flow rate, bottom flow 
rate, stopsiproduct ‘low rates neds xeele OWemrda tc SUC aim LOW 
rate, condensate level, top product composition and bottom 
DrOGUCtaMCcCOMmMpOSit10N, Were, recorded ~on Sdisc. using DISCO. 
These variables plus the column pressure were monitored 


every 3 minutes by a safety monitor program. For implement- 
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ation of the control scheme, only the manipulated variables, 
Steam and reflux flow rates, were under a supervisory type 
of control¥with the outputs*of “the control calculations be-= 
ing used as the setpoints for the local analog PI control- 
lers while the other variables were under local analog PI 


cont rod’ 


6.3 Implementation of Control Algorithms using DISCO 
Bhefoverarivone afoDvisc0seisia@ibasied won aAheh concept fof 
table-driven processing. One of the distinct features of 
DISCO 1istitsr7T1exibletdata pbaset structure. thachtlactivityafor 
loop) of DISCO is buiit by adding different segments which 
perform different tasks, to the activity header which iden- 
tifies itself from other activities. The operation within an 
activity is task-driven, i.e. different segment processors 
(or FORTRAN subroutines) will be called by the DISCO main 
program to process the data base of a particular activity. 
For a single input-single output feedback control loop, the 
Simplest structure of an activity contains an activity head- 
er; an input segment which acquires the measurement of the 
controlled variable; a control segment which performs the 
control calculation according to the control algorithm cho- 
sen; and an output segment which sends the calculated value 
of the manipulated variable to the process. A description of 
the organization and operation of DISCO is given by Brennek 


(5) and in documentation available from the DACS Center (62, 
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685.64 soo nae PhLiSstseotaethe sactivities:gusedednfthishwonk lis 
available from the DACS Centre. 

During the project, the existing input and output seg- 
ment processors (5) were rewritten to allow for greater 
flexibility for the users. A control segment processor which 
eongepe nhommscne gcongye be calculation using Equatione 3.161. 
was written to allow for the implementation of different 
control algorithms using the same activity. A special seg- 
ment processor, the program scheduler, was also written for 
actuating the gas chromatograph operation or data acquisi- 
tion. The program coding and the user's manuals for these 
segment processors are available from the DACS Center (62, 
63,64 65). 

A set of interactive programs which can be used to di- 
rectly access the data base of an activity were written to 
provide for implementation of on-line control algorithm tun- 


ing. These programs are available from the DACS Centre. 
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Since the, distal lation scolumiasis. af highl yeinonilsinear 
system, especialiy the response of bottom product composi- 
tione tos-step changes.im feed flowsrate, iGbs1S necessary for 
a control algorithm to use different sets of tuning con- 
stantsiefor controihingn datiierent edusturbances suntordenito 
yield optimal performance. In the experimental evaluation, 
the performance of the different algorithms for a -25% step 


changeadintfeede flowfratencis of particularm!interest., since 
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this change is the most severe disturbance for the pilot 
Scale distillation coiumn. Summaries of the performance of 
the different control algorithms, namely the digital PI and 
PiDgcontrol valgorithm, the Dahlin wand Smith tpredictor con- 
Geol ‘algorithms Aavyareg givenoin Table 6.4.1 wand 604s. 2% athe 
bottom composition control responses using the above algo- 


mucimseare shown in higures 6.4.1 to 6.4.8. 


655 Discussion of Regulatory Control Results 

Using the IAE value as the performance criterion, it 
was found from the experimental evaluation of the different 
control algorithms=tor a -25% step change in fe.d flow rate 
that the Smith predictor provided the best performance as 
indicated by its lowest IAE values (cf. Table 6.4.1). How- 
ever the controlled behaviour achieved using the Dahlin 
algorithm was virtually equivalent as indicated by its IAE 
Vahuest(cf. fable 6.4.1). Tests using the digi,talisPID Faligo-= 
rithm showed that by careful tuning, its performance was 
comparable to that achieved using the Smith predictor and 
the Dahlin algorithms. The performance of the digital PI 
akgorithm was infervor to that. of atnesocher al goniunmsnerds 
indicated by its highest TAR values (cf. Table 6.4. i). 

Although the Dahlin algorithm is not specifically a 
time delay compensation control scheme such as the Smith 
predictor scheme, which cancels the time delay term in the 
System characteristic equation, its successful performance 


compared to the performance of the digital PI algorithm 
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Tab Vero..47 1 


Summary of Experimental Results for 
Contvolmor Bottom Composition tor 
a -25% Step Change in Feed Flow Rate 





Algorithm So Ky Kp Kp Tp Tq IN Figure 
PI -0.607 -0.095 — — — — 5148 6 4.1 
BID y -0.523 -0.117 -0.200 — — — 4477 6.4.2 
Dahlin \=164s — — =o 4 9167187542535 6" 4.53 
Si imnmC Cee Ol A) ei U2 toe 2 7088/2 18714236) 6.40.4 
Ke = (g/s)/wt.% Kp = wt.%/(g/s) 
Kr = (g9/s)/(wt.%-s) Tp = S 
Ky Sg (Si) /' (Wie oS ) lees 
IAE = wt.%-s 

Table 6.4.2 


Summary of Experimental Results for 
Gontrol=of Bottom Composition for 
a +25% Step Change in Feed Flow Rate 


eR ES 


Algorithm K 





Pal -0.551 -0.122 — — — — 4014 6.4.5 
PID -0.503 -0.125 -0.075 — — — 4158 6.4.6 
Dahlin a=231S — — =o [4 9185 Ley 502286. Ace 
Smith Pred. -0.711 -0.185 — =F 0S / 2 Luh 06 2a 6 406 


Ke = (g/s)/wt.% Kp = wt.%/(g/s) 
KT = (g/s)/(wt-%-s) Ty = § 

Kp = (gs) /iwt-%/s ) ees 

IAE = wt.%-s 
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results from the consideration of the process characteris- 
tics during the controller design procedure. 

For the PID control algorithm, properly tuned deriva- 
tive action definitely provided significant compensation for 
the time delay as indicated by the much lower IAE values 
fofi2 babke 26 44: l)sethana thoseainesuiliting fromercontrolli suis iing 
the PI algorithm. 

iInvterre Icasceeof Pbottom composi tivonecontr oll 9for trans+25% 
step change in feed flow rate, the control performance using 
the Smith predictor, Dahlin and PID algorithms was inferior 
to that achieved using the PI algorithm as indicated by the 
IAE values shown in Table 6.4.2. Some explanation for this 
performance is the fact that these algorithms were not tuned 
for this particular disturbance. A further contributing fac- 
tor to this performance is the fact that when the feed flow 
rate is increased, the distillation column operates close to 
the flooding region so the assumed transfer function models 
for the distillation column dynamics are probably not reli- 
able. Consequently, the control pertormance for an increase 
in feed flow rate do not provide a valid evaluation of these 
algorithms. 


6.6 Performance of the Algorithms for Servo Control 


Since load disturbances occur more frequently than set 
point changes, most of the controllers in the process indus- 
tries are tuned for regulatory control. From this point of 


view, it was more practical to investigate the performance 


eo] > 
ai an a 2369 ail 5 ietde%ebi ete gfj mort etfu2o7 | 
seb velfoxndnoa sat paryyo zor? ; 
miirvopis foaxtnes G19 snd 304 _- A 
’ 5 ryorg eletintteb nors96 syitt | 
\ 


yaloab omf3 saz 


- 
. 
~~ 
w& 


s20n: ana (f.e.a ~ofdst tan s 
mitinoptss149 ahs 
=a to) «s26o sfF aly a 
6a wolt bost nf eensnd qat2 
—~ ia hn bins ywJIstestg Ajime ght 
ry f G ont f2 bsvainos tend ot 
omni SS. ave & 5 yi nwofre 2euTsy WAI * 
zoo Jans 3 gn3 27 9]nsm10t78q” 
u 
‘ TV UT “¢ anwed vy : ‘6 uorgI76q 2rns 10T | 
| — 
ise? otf sowemrottieg 2tdd of Tem 
ir . ,bazse7on! ef 3361 | 
7 ~ ys 
\a bene suuzet ast? de notaet entboolkt ant 
tien ton utd 7 saulos wottafitaetb, add 7aRe 


S26s19hi HS 101 namrorteg Lownie. si ‘UT Smaupsengy .sid6 a 
4 - =e 


ated TO AOTMbulsvs afhisy 6 sdivowg Fon OO Sie" woft bay? nF 


= | .emastropii 





120 


Or ethers avi ferent controliaaigonithms for servo control using 
the tuning constants?’ Used? forttreguiatonry control than to 
Gunesthe algorietnmSestrictiveroe Sservorcontrol operation. In 
this work, the tuning constants used for composition control 
with the column subjected to a -25% step change in feed flow 
rate were also used in the algorithms for controlling the 
bottom composition for +1% changes in set point. The IAE 
values for the performance of the different algorithms are 
Seven ne suabhkese 605.1 and 6.5.2. she composition responses 
using the different control algorithms are given in Figures 


Opp-1 £0,6-5.8. 


—_— 
—_— ee eee, eee} 


From the results shown® ini elables 6.5.1 and 6.5.2, it 
can be seen that the Dahlin algorithm with the same tuning 
Canstants as usedisiormerequiatoryi control’ "ior a -25% step 
Change in feed flow rate, provided the best and the most 
Gonscustent performance for =both +1% step changes inset 
point. The overall performance of the Smith predictor was 
Ghose Coecthati.of ythel Dahil in—al gorithm. spank cuvvarly)/ ion Sthe 
-1% step change in set point which gave an IAE value nearly 
equal to that achieved using the Dahlin algorithm. Control 
performance using the PI and PID algorithms with constants 
determined for a -25% step disturbance in feed flow rate was 
inferior to that achieved Susing jtnhe” Dahiin, algorithm ‘or 
Smithepredvccor.. |NiS = tsesinaicaved=not only by. thes higher 


IAE value (cf. Tables 6.5.1 and 6.5 2), but by the composi- 
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THONMEGeSpONSesepsShowtieswn Figuresn Gia5 Af 6r 5628856 25125 Tand 
6.5.6 which are more sluggish and oscillatory than those 
that resulted when using either the Dahlin or Smith predic- 
SOnceu gonpahms etc Teuraigunes, 67563 .dec:.oe4?, 62 567 danidye6). 5t She 
Furthermore control using the PI algorithm for a +1% step 
enangerin setaipoi ntealof ~ehigune te5 oS horesulted pin suchran 
oscillatory response of the composition that such behaviour 
Cannot considered to be acceptable. The IAE values in Tabies 
O.5-1 Gdnd's61.ibs 2enshowede tivati the tPikDi algorathms with, ats 
derivative action provided a significant improvement in con- 
trol performance compared to the PI algorithm for control- 


ling processes containing time delay. 


6.8 On-Line Tuning Techniques 

bheachoicehior wheiginitiala set ofscontrolleruconstants 
TsHimpontant for anyutype of control sappkication; withsthe 
selection of the initial tuning constants dependent on the 
user's knowledge of the process and the tuning technique. 
From the simulation studies, employing the improved tuning 
technique as outlined in Chapter 5, forthe feed flow rate 
disturbances, it was found that the values of the model 
parameters used in the Dahlin algorithm were approximately 
the average values of the model parameters of the load and 
process transfer functions. 

However, in the case of the online tuning of the Dahlin 
algorithm, the initial values of Ky and the Tq used in the 


algorithm were chosen to be the same values as those found 
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in the open loop transfer function of the process; the ini- 
tial value of Tp in the algorithm was set to approximately 
the average value of the time constants in the load and 
process transfer functions. The desired closed loop time 
constant A was initially set to be approximately 80% of the 
Tp of the process model used in the algorithm. From experi- 
ence, these parameters ensure a stable start-up of the con- 
trol loop. The tuning of the Dahlin algorithm proceeded by 
first adjusting the value of A and then the values of Ky and 
ae A summary of the results obtained using the Dahlin algo- 
rithm is given in Table 6.6.1 and the composition responses 
amcagiven? in Rhitgunes_6%6 -teto <6 76,29). 

As shown by the IAE values in Table 6.6.1, the perform- 
ance of the Dahlin algorithm improvea as the value of i was 
decreased. Furth.r minimization of the IAE value of the 
Dahlin algorithm resulted by applying the improved tuning 
procedure. As shown by the IAE values in the Table 6.6.1, 
the control performance improved gradually as the value of 
Ss was decreased while the value of Tp was increased. From 
Brqures o676e5e tonet 649 eand) 6.46 3.0m tiheaniaberseentihat lone 
initial magnitude of change in the manipulated variable, in 
reaction to detection of the error, increased as the value 


of K) was decreased. This behaviour is explained by the fact 


that the Dahlin algorithm, by virtue of the low value of Kp 
in the transfer function model (meaning that the bottom com- 
position is insensitive to changes in steam flow rate) cal- 


culates a large correcting action in order to force the pro- 
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Table 6.6.1 


Summary of Experimental Results for Dahlin 
Algorithm Control of Bottom Composition, 
for a -25% Step Change in Feed Flow Rate 
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cess to follow the desired closed loop response. Further- 
more, it can be seen that the sluggishness of the composi- 
tion response disappeared as the value of Ty in the Dahlin 
algorithm was increased. The reason for increasing the value 
of Ty is analogous to the strategy of decreasing the value 
of KD. By virtue of the large value of Tp meaning that the 
Dahlin algorithm is controlling a slow process, the frequen- 
Cy of change in the manipulatec variable has to be higher in 
order to force the process to follow the desired closed loop 
response as dictated by the alyorithm. 

In -the ‘case of “Cuning. them sm icheprediclOm mat heusiunad 
set of model parameters used in the Dahlin algorithm were 
also used in the Smith predictor therefore the tuning of the 
Smith predictor invoived primarily adjusting the value of Ko 


and KT. A summary of the experimental results achieved with 
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Figure 6.6.1 Dahlin Algorithm Control of Bottom 
Composition for -25% Step Change in 
Feed Flow Rate (Run E-DAHO1; 
Ko = -5.63, Tp Sao Ue te 1 87. 0), 
A= eO JU) ) 





13 


1B. 68 13.64 16.648 


18. 84 


Pha resetrnet | Ronn sre pores a | 


=> 


fo 


1. HS 


ree eP at 4 
be, 
Sh 
a, 
Fe. 
H i 


ETT T IT 
nS LET ST rf 


Figure 6.6.2 


134 


1 
a. yee tn, aries Tete neta he Mie ee oe 
= PY egy ret 

32 


| | 
$6. 88 = oe. BE 1h, Ge 
TIRE “SEE Cale 


7aE | 
vr, - Wey ter 
1 BG Ae ae aS 
m ATED Tare Ti AE 





Dahlin Algorithm Control of Bottom 
Composition for -25% Step Change in 
Feed Flow Rate (Run E-DAHO2; 
Ky = -5.63, Tp = (O00, iq = 187.0, 
A" =. 863.0) 
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Dahlin Algorithm Control of Bottom 
Composition for -25% Step Change in 
Feed Flow Rate (Run E-DAHO4; 
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Figure 6.6.5 Dahlin Algorithm Control of Bottom 
Composition for -25% Step Change in 
Feed Flow Rate (Run E-DAHO8; 
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Dahlin Algorithm Control of Bottom 
Composition for -25% Step Change in 
Feed Flow Rate (Run E-DAHO9; 


K, = -4.86), Tp > 289.0, Tq = 187.0, 


¢° 





ue meee “ - 8 


‘, 
4 
ee 


CWT 


E 


¥) 
ee 


39 


12,84 


£ 

= Pea, a arhe rahe ime ay ae _* ‘ 

4, us ore ral bat Be Fe Sema ht ‘ fee ane! 4a! area ar ates encase 
aoe 
1 


18, Ag 





Ge 


he 


~. -, 
. TiarTY f aS Pome me Sa Trottier. of oa rome 
aa iy hy pg ol 7 ti presse oti Te See aie AT Tyr 





Figure 6.6.7 Dahlin Algorithm Control of Bottom 
Composition for -25% Step Change in 
Feed Flow Rate (Run E-DAH10; 
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Summary of Experimental Results for Smith 
prediretom Contrall of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 





Run Ky Tp Pa Ke Ky IAE Figure 
E-SLP09 3 sf496918 187 9050.8300 209209 446,09 6 6.10 
E-SLP10 =301467918 187 3020.869.9 205209 44361 G65 12 
E-SLP19 =3eR4BE918 1829150. 71L0.400185 44549 §-G.12 
E-SLP23 =O 52, 98072" 1500870.74209204213 442316 6.4.4 
Ke = (g/s)/wt.% Kp = wt.%/(g/s) 

Kr = (g/s)/(wt.%-s) Tp = 8 

IAE = wt %-s Woe Ss 


composition control responses are given in Figures 6.6.10 to 
Gao. 12" 

The tuning of the Pl@and PID algorithms was basically a 
trial and error procedure starting with the initial set of 
tuning constants. calculatedh usangiythesecornehatians of Smith 
(51). Summaries of some of the experimental results are 
Gavenm sin ‘Tables 6.6.3 and 1.6.6.4. Ihe bottom composmtiion 
responses obtained operating under PI and PID control are 
snowmein Figures 6646.13 Co 6.6.20 

In common with any tuning procedure, the improved tun- 
ing procedure is a trial and error routine so the amount of 
time required to obtain ~aesatisfactocyeset sotsstuninge con - 
stants is dependent on the process being controlled as well 
as the experience of the individual performing the tuning. 
During the experimental evaluation of the tuning procedure 
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Summary of Experimental Results for PID 
Algorithm Control of Bottom Composition 
for a -25% Step Change in Feed Flow Rate 





Run Ko K] Kp AE Figure 
E-PIDO5 -0.636 -0.200 -0 666 5236 a cllsye: KS) 
E-PIDO9 -0.701 -0.105 =0.333 4554 6.6.14 
E-PID13 -0.573 -0.100 -0.200 4558 GrlO tO 
esp IDT5 -0.659 -0.110 -0.170 4986 Eig tes dls 
E-PID21 -0.523 = (Quel 7 -0.200 4477 6.4.2 


pret (o/s) /wons 
T=) (g/s)/(wt.%-s) 
(g/s)/(wt.%/s) 


=] 
iT} 


Table 6.6.4 


Summary of Experimental Results for PI 
Algomitnme Gontrol of Botcom Composition 
0G sar -copeOLepeLiange ineheed@Elow Rate 





Run Ke K] IAE Figure 
E-PI14 = Or Ona -0.122 5850 OpeOrdy/, 
E-PI18 -0.559 S/N Saba 5364 62 O.u6 
E-P1I20 =0 55a -0.106 6534 OO 
E=Pi22 -0.578 -0.100 6282 6. Gine0 
E-P124 -0.607 -0.095 5148 6.4.1 
Ke = (g/s)/wt.% 

Kies 
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IAE = wt.%-s 
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Figure 6.6.13 PID Algorithm Control of Bottom 
Composition for a -25% Step Change 
im feed. Flow “Rate “Run “E-PIDOS 
Kc= -0.636, Ky= -0.200, Kp= -0.666) 





OTS, 
“SD 


q Lae’ 


ST 


Ae 
a 4 


OWT: 


Et 


=, 


16. 64 


18, ag 


Aj 


i 


1B. 


STS) 


148 


tyr ss ; 
4a Wy 4 
L ? m 
wok we 
i arch anh 3 ea ae all 4 
| P H ihe hE Veet Parke Eh mal 





rane 
| 
| 
| 
} mee try TM = 
Ty mytatiqer, al Cys Tr i, mT, sit CAT Tj 
Wine! oe 1 9 -5T) mm a th pe saan i US “ae 
- = - mT 
iT} or 
3 a all) 
ral) a nae 
a) 
m a 
T z 
t } J RE, a sah re wigs ae 
a. oe oe, A 46. 08 68. 6 St cae 166.68 _ 
ievaies teeta) CHLbee= 


Figure 6.6.14 PID Algorithm Control of Bottom 
Composition for a -25% Step Change 
in) Feeds -Fitow, Ratesr (Run E-pIp09. 


Keay OR PONY, Ky=, -05105, Kp=. -05 333) 





ar 


CWT 


oy 


1B. f4 


fi 


i ft 
ey A Sin ‘ 


at we ‘y, wo anaes ea AAD ne MA atte 


o 
| ITI 17} —/lhtar, v 
iTiret- ace in minal rrert te Ture. THT eer, 
Che Ty ATT Ee, = Ti Ts Tete tyre te har Ta 
. “GH ' ae ho iS eae & 
aor 


a a 





TIME KSECD 
Figure 6.6.15 PID Algorithm Control of Bottom 


Composition for a -25% Step Change 
in feed flow Rate’ (Run '£-Pibis: 


Kee -0.573, Ky= -0.100, Kp= -0.200) 


149 





7. 


deen a7 Gent 


i 


ah 





Pome 
Las = e | 


ar ¢ 
AY 


CTS 


16.68 


13 
at, 


1B. Ay 


=, 
= 


Zar 


150 


- xi 4} 
4 7 its, ee . ‘ 
pce beng L— a hey al, age aa Pe Jong 


Se cat \ 


t [ hy ae as ae! - pao al a 





a 
Le Mm mete, 
~ a es rerfitTeTs uy 
AD ity 7; iT . Giri td 4 teres rears AOTin, 
ud i 
rar} Oe. Tit 
7 Tit 
a Th 
tT m7 


; 


a 


Figure 6:6:16 PID “Algorithm Control of Bottom 
Composition for a -25% Step Change 
in Feed Flow’ Rate®—-(Rinm E=Prols; 


Kos -0.659, Kyr= 20310 Kp= -0.170) 







Ss == aa 


7 7 TEEN or, 
PR trea Ps. eat 2 eo 





aS) 


ST ¢l 


. BH 


16.64 


is 


15a 


18.88 


, 
s 


ry be 


aT 


18. 


. 68 


rd 
t 


AES 


4. 88 





ages 
| h ent 
1! ro lee r. eon, Fs 4 r eis 
‘ { d % » oan hee = a abt eee Feyenaar Meee a ¥ 
L A j : oe . htt tr 
i) | 
: 4 
a Mond 
eee 
BB 28. Ba 4h. aa Bi. Ge SA. We L6H. GF 
Tine @seco CAB? = 
Grr 
{ iti 
Thr, 
=. iT 7) om in 
pots, ne Mig: Cree aytern Tn 
is “ m1 { Inyti— SS SS 
m % ~ mi 
| ni 





Hoe 2G. AG S68. Ge ca, Be SS Lee. be 
TIME (SEC? cmp = 


Figure 6.6.17 PI Algorithm Control of Bottom 
Composition for a -25% Step Change 
in Feec Flow Rete (Run: E-P iia: 
Ke = -0.551, Kr = -0.122) 
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amount of time required to tune the Dahlin algorithm was 
comparable to the time required for the PID algorithm to 
yield a similar control performance. However, the Dahlin 
algorithm would tend to be favoured by control engineers 
because of the nonoscillatory responses of both manipulated 
and controlled variables (cf. Figure 6.4.3) as compared to 
those that result from using the PID algorithm (cf. Figure 
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7. Conclusions and Recommendations 


From the digital simulation studies, the following con- 


clusions can be drawn: 


ie 


The nonlinear distillation column model as developed 
by Bilec (4) and subsequently moditied in this study 
Provided a satisfactory representation of the dynam- 
ics of the distillation column. 

First order time delay transfer function models pro- 
vide a satistactory representation of the column dy- 
namics for control purposes. 

Pnetne case=of set point changes, the deadbeat; Dah- 
Tin and Smith predictor algorithms with only minimal 
model errors provided improved control performance 
compared to that achieved using the discrete PI and 
PID algorithms. 

For changes in feed flow rate, especially when the 
dynamic of this disturbance is much faster than that 
of the process, the deadbeat, Dahlin and Smith pre- 
dictor algorithms tuned by using conventional tuning 
methods failed to provide satistactory performance 
compared to that achieved using PI and PID 
algorithms. 

Alithough the deadbeat algorithm does not have any 
provision for on-line tuning of parameters in the 


conventional sense, the process model parameters 
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used in the algorithm can be used as on-line tuning 
parameters. 

The above strategy can also be applied to the Dahlin 
and Smith predictor alyorithms especially when these 
algorithms are applied for regulatory control opera- 
Ciens? 

US em oTmedenivati vemactiion, Pethatwe isatusing:! the? tpn 
algorithm instead of the PI algorithm, significantly 
improved the control performance for either servo or 


regulatory control operation. 


From the experimental results, the following conclu- 


Sions can be drawn: 


ha 


Fon | ynequilatoryss controdst eitheyncontrol performance 
achieved using the Dahlin and Smith predictor algo- 
rithms tuned by using the improved tuning procedure, 
was much better than was possible using the PI algo- 
rithm. However, with the proper tuning, the pertorm- 
ance achieved using the PID algorithm was found to 
be comparable to that achieved uSing the Dahlin and 
Smith predictor alyorithms. 

For set point changes, control using the Dahlin and 
Smith predictor algorithms with the tuning constants 


used for regulatory control provided superior con- 


‘trol performance to that achieved using the PI and 


PiDealbgor itihm s 
The results confirmed the simulated results that the 


control behaviour using the PID algorithm is much 
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better than is possible using the PI algorithm for 
both servo and regulatory control. This demonstrates 
the effectiveness of derivative action on the con- 


trol of processes containing time delay. 


Recommendations for future work suggested by this study 


include 


ne 


the following: 

From the design procedure for the Dahlin algorithm, 
it can.be seen that this algorithm is designed pri- 
Mapihye for, servo canenol cope ga tion.) A litthough) dite proz 
vides satisfactory performance for regulatory con- 
Cro lmeopeiat FOnmeINNnenIe@tte TSiatunedsgeby fusing ethetam= 
proved tuning technique, it is more appropriate to 
include@ thes loadl*transfer 'tunction tmodele imi’ the 
design procedure if regulatory control operation is 
the primary concern. 

oamid amy, a itgatheiesmith predictorsts used’ fore con= 
trolling systems subjected to frequent load distur- 
bances, load prediction as suggested by Meyer (38) 
Should be incorporated in the control algorithm. 
Other digital control algorithms which employ design 
techniques that are similar to the direct synthesis 
method and have thes formeot sequata cies sit sucimards 
the Kalman alyorithm (10), should be studied and 
compared with the-PiD aluorithm=to access their per- 
formance for controlling time delay processes. 

Vences tie deadpedu sl aiimn candesmitiepredictor con - 


trol algorithms are model dependent, the control 
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performance of these algorithms using a higher order 
Cyanster mm UnClLON= Modell In sthesa! gorithm synthesis 
Should be studied. 

Since the performance of these alyorithms depends on 
the accuracy of the parameters of the process trans- 
fer Tunction modely it as*suggested’ that® an on-line 
parameter estimation scheme which will operate in a 
Supervisory mode should be incorporated into the 
algorithms. This scheme could easily be implemented 
using DISCO. At each sampling time, the parameter 
estimation scheme would calculate the parameters 
used in the transfer function model and the coeffi- 
Cients of the control algorithm will then be calcu- 
lated subsequently and transferred back to the DISCO 
data base. Then, the control output would be calcu- 
lated based on the new coe,ficients and the error 
and output sequences. The advantage of this approach 
compared to that uSing a parameter invariant model 
iss uthat “itwswould account. for changes in process 


dynamics due to disturbances. 
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